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SOMMAIRE 


Des  profiles  de  vitesse  et  de  la  fluctuation  des  vitesses  en 
regime  turbulent  ont  ete  obtenus  experimental lement  dans  la  couche 
limite  et  la  region  central  de  tuyaux  circulaires  a  l'aide  d'une 
methode  photographique  utilisee  pour  mesurer  les  vitesses  i nstantanees . 

Les  donnees  experimentales  ont  permis  de  deceler  des 
erreurs  d ' i nterpretation  des  profiles  de  vitesse  pour  des  systemes 
dans  lesquels  une  reduction  de  perte  de  charge  est  observee,  et 
ont  offert  un  appuie  au  phenomene  d' ejections  comme  mechanisme 
responsable  de  la  reduction  de  la  perte  de  charge. 

En  particulier,  il  a  ete  demontre  que  la  correlation 
empirique  presentee  par  Bogue  doit  etre  incorporee  au  profil  semi- 
logari thmique  de  vitesse  de  Seyer-Metzner  pour  que  les  donnees 
experimentales  soient  representees  par  des  courbes  parallel les  a 
I'equation  de  Prandlt. 

II  a  aussi  ete  demontre  que  la  fonction  B,  definie  par 

Seyer  and  Metzner  et  presentee  a  la  figure  11-12,  peut  etre  utilisee 

+ 

dans  les  limites  de  faible  valeurs  de  e  ,  mais  que  la  valeur 
assymptotique  de  32  n'est  pas  unique. 

II  est  apparu  evident  que  la  fluctuation  de  vitesses  des 
systemes  avec  diminution  de  perte  de  charge  peut  etre  plus  grand  ou 
plus  petite  aux  systemes  Newtoniens  selon  la  base  de  comparaison 
choisie.  Les  donnees  experimentales  ont  demontre  que  le  parametre 
d'echele,  definit  par  Rudd,  n'est  pas  aussi  general  que  desire. 
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ABSTRACT 


An  experimental  study,  based  on  streak  photograph-deter¬ 
mination  of  instantaneous  velocities,  was  directed  at  determining 
the  velocity  and  intensity  profiles  within  the  boundary  layer  and 
core  regions  of  circular  pipes. 

The  measurements  resolve  several  discrepancies  in  inter¬ 
pretation  of  earlier  velocity  measurements  in  drag  reducing  systems 
by  lending  support  to  the  ejection  phenomena  as  the  mechanism 
controlling  drag  reduction.  In  particular,  it  is  shown  that  Bogue's 
empirical  correction  function  must  be  incorporated  in  the  semi- 
logarithmic  velocity  profiles  to  quantitatively  represent  the  data 
of  drag  reducing  systems.  Failure  to  do  so  has  led  to  confusion  in 
determining  the  proper  form  of  similarity  law. 

Study  of  the  details  of  Seyer  and  Metzner  velocity  profile 
intercept  function,  defined  in  equation  11-12,  confirm  the  function 
for  e+  up  to  approximately  12,  but  show  the  estimated  maximum  value 
of  32  is  not  unique. 

Turbulence  intensities  may  be  viewed  as  being  higher  or 
lower  than  Newtonian  intensities  depending  on  the  basis  of  comparison. 

The  data  show  that  Rudd's  scaling  factor,  which  superimposes  his 

*  ★ 

data  for  fixed  U  ,  does  not  account  for  changes  in  U  and  is  there¬ 
fore  not  as  general  as  earlier  anticipated. 
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INTRODUCTION 


Considerable  numbers  of  theoretical  and  experimental 
studies,  in  the  last  10  years,  have  been  focused  on  the  mechanism  of 
turbulence  in  an  attempt  to  reveal  the  character  at  the  wall 
responsible  for  drag  reduction  observed  in  the  presence  of  polymeric 
additives.  In  principle  a  detailed  knowledge  of  the  flow  in  the 
wall  region  would  aid  one  in  choosing  the  correct  constitutive 
equation  or  variables  to  correlate  drag  reduction.  Many  mechanisms, 
as  well  as  correlations  of  friction  factor  and  velocity  profiles, 
have  been  proposed.  These  correlations  represent  the  available 
data,  but  none  have  been  tested  very  severely  over  a  wide  range  of 
variables.  In  particular  the  available  velocity  measurements  in 
drag  reducers  do  not  provide  a  detailed  verification  of  the  semi- 
logarithmic  velocity  profile.  Also  great  confusion  is  found  from 
a  survey  of  the  measured  relative  intensities  owing  to  the 
uncertainty  in  accuracy  of  the  results  and  in  how  to  scale  or 
compare  them  for  different  systems. 
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1-1  MECHANISM  OF  TURBULENCE 

Since  Prandtl  developed  the  'Boundary  Layer  Hypothesis'  in 
1904,  many  studies  of  properties  of  Newtonian  turbulence  (19,21,14, 
28,45,16,18,38)  have  indicated  that  the  character  of  the  flow  in  the 
wall  region  was  responsible  for  most  of  the  creation  and  dissipation 
of  the  turbulent  energy. 

It  is  only  recently  that  a  detailed  physical  picture  of  the 
mechanism  was  presented,  mainly  by  the  work  of  Bakewell  and  Lumley 
(3),  Kline  et  al  (20 )  and  Corino  and  Brodkey  (?).  Bakewell  and 
Lumley  indicated  that  the  dominant  large  scale  structure  of  the  flow 
in  the  boundary  layer  consists  of  randomly  distributed  counter¬ 
rotating  longitudinal  eddy  pairs  of  elongated  extent  in  the  flow 
direction.  The  structure  of  these  eddies  was  inferred  from  space- 
tirne^ correl ation  functions  of  the  fluctuating  velocities.  The 
streamlines  of  these  eddies  were  pictured  as  to  push  low  momentum 
boundary  layer  fluid  toward  the  core  region,  so  that  the  renewal 
of  fluid  was  by  flow  in  the  circumferential  direction. 

Qualitatively  similar  patterns  of  flow  in  the  wall  region 
were  visualized  by  Kline  et  al  using  a  dye  ejection  technique  and 
in  more  detail  from  tracer  photographs  by  Corino  and  Brodkey. 
Ejections  of  fluid,  originating  from  a  low  velocity  region  adjacent 
to  the  viscous  sub-layer,  were  responsible  for  extracting  energy 
from  lumps  of  fluid  originating  in  the  main  flow  and  converting  it 
into  turbulent  energy.  They  observed  that  these  ejections  were  of 
large  scale  and  moved  through  the  boundary  layer  until  broken  down 
to  small  scale  turbulence  by  mixing  with  the  main  flow. 
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The  simple  mechanism  of  turbulence  then  is  that  the  turbulent 
shear  stress  is  generated  by  the  radial  transport  of  low  momentum  fluid 
by  the  large  eddies  or  'bursts'.  The  magnitude  of  the  stress  is 
determined  by  the  rate  of  fluid  transported  radially  as  well  as  its 
axial  momentum.  The  rate  of  fluid  transported  in  turn  depends  on  the 
frequency  of  occurrence  of  the  eddies  as  well  as  their  size.  The 
visual  studies  of  wall  turbulence  provide  a  starting  point  for  inter¬ 
preting  turbulence  measurements  in  drag  reduction. 

Many  mechanisms  have  been  suggested  in  order  to  explain  the 
drag  reduction  for  flow  of  viscoelastic  fluids  in  pipes  (32,39). 
Generally  it  is  agreed  that  elasticity  of  the  fluid  is  responsible  for 
the  drag  reduction.  Furthermore,  as  emphasized  by  Wells  and  Spangler 
(48)  who  showed  by  injection  of  polymer  in  the  boundary  layer  that 
drag  reduction  is  controlled  by  the  flow  in  the  wall 'region,  the 
elasticity  must  modify  turbulence  in  the  wall  region. 

From  a  continuum  point  of  view,  Seyer  and  Metzner  (17)  have 
shown  that  the  drag  reduction  can  result  from  the  resistance  to  flow 
caused  by  stretching  in  the  large  eddies  or  'bursts'.  Although  their 
arguments  are  based  on  a  steady  state  analysis  of  eddies  as  pictured 
by  Bakewell,  Denn  and  Porteous  (S  )  have  shown  identical  ordering 
parameters  appear  if  a  transient  analysis  of  the  eddies  is  carried 
out.  It  is  difficult  to  say  what  the  effect  of  resistance  to  flow  is 
on  the  structure  of  the  eddy.  For  example,  the  eddies  may  be 
increased  in  size  which  would  appear  as  a  thickening  of  the  boundary 
layer  or  they  might  be  reduced  in  frequency.  These  points  of  view 
also  have  resulted  in  dimensionless  groups  which  are  the  same  as  for 
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stretching  arguments  (24,17)  but  do  not  show  conclusively  what  altera¬ 
tions  in  the  structure  have  occurred. 
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1 1-2  SIMILARITY  LAWS  AND  MEASUREMENT  METHODS 

For  drag  reducers,  apart  from  completely  empirical  formula¬ 
tions,  most  studies  have  sought  to  modify  the  Newtonian  similarity 
laws,  so  that  they  described  drag  reduction.  For  Newtonian  fluids, 
Millikan  (27),  for  example,  has  shown  from  a  dimensional  analysis 
that  the  velocity  profile  law  is 

U+  =  A  In  Y+  +  B  1-1 


where 


+  _  * 

U  -  yu 


yU*/v 

yp 


1-2 


t  ,  =  DAP/4L 

w 

The  generally  accepted  values  for  A  and  B  for  Newtonian  fluids  are 
( 5  )  A  =  2.46  and  B  =  5.6 

Equation  1-1  is  an  accurate  representation  of  all  consistent 
Newtonian  data  for  Y+  less  than  400,  but  considerable  deviation  occurs 
for  larger  Y+.  Correction  factors  have  been  empirically  found  (13, 
5,2 7).  Bogue's  correction  factor  represents  the  best  fit  of  avail¬ 
able  data  and  is  defined  as 

.  U-0-8)2 

C3U.f)  =  0.05/|  e  °'15  1-3 


where  £  =  y/R 
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and  the  friction  factor 


t  /I  —  2  T  . 

f  =  \'z p  'V 

With  Bogue's  correction  factor  the  form  of  the  velocity  profile  for 
Newtonian  fluids  is  then 


U+  =  a  In  Y+  +  B  +  C3U,f) 


1-5 


From  Equations  1-4  and  1-5,  the  friction  factor  relation  for  Newtonian 
fluids  follows  directly 


1 

/f 


=  A  In  (NRe/f)  -  0.4 


1-6 


For  fluids  represented  by  a  power  law  model  in  which 


8<V  n 

T«  ■  K' 


1-7 


Dodge  and  Metzner  (10)  formulated  the  friction  factor  correlation 
based  on  the  generalized  Reynolds  number 


1 

/f 


^05  l09  <NRe  f 


1- 


n 


)  - 


(n‘) 
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where  the  generalized  Reynolds  number 


N 


Dn ' <u  >2-n* 
Re =  — 


K  8 


Notice  Equation  1-8  is  simply  a  generalization  of  the  Newtonian 
Equation  1-6  to  account  for  the  non-Newtonian  properties.  Similarly, 
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Bogue  (5  )  measured  velocity  profiles  for  pseudoplastic  fluids  and 
formulated  the  velocity  profile  equation  such  that 


U+  -  C3U,f)  =  5.57  log  (Y+)1/n'  +  I(n',N^) 


1-70 


where  I  is  a  weak  function  of  n',  Nn  . 

Equation  I - 1 0  when  plotted  gives  a  family  of  parallel  curves 
to  Newtonian  fluid  with  slope  A  =  2.46  for  the  core  region  with 
intercept  higher  as  n1  decreases  and  N^g  decreases.  He  also  found 
that  turbulent  velocity  profiles  for  purely  viscous  non-Newtonian 
fluids  are  essentially  the  same  as  those  for  Newtonian  fluids  when 
normalized  to  the  mean  velocity. 

For  drag  reducing  fluids  it  is  obvious  that  one  of  the  major 
effects  appears  as  an  increased  intercept.  Based  on  Pitot  tube  measure¬ 
ments,  Elata  et  al  (11)  modified  the  law  of  the  wall  by  i ncorporati ng 
an  extra  logarithmic  term  to  interpret  the  higher  intercept  such  that 


U+=AlnY++5.5+al  In  [ 


1-11 


In  this  equation,  A  is  assumed  the  same  as  for  Newtonian  fluids  and 
the  effect  of  the  elasticity  of  a  solution  is  expressed  by  the  relax¬ 
ation  time  e  which  is  assumed  constant.  The  parameter  a1  is  an 
empirical  function  of  the  polymer  characteristics  and  the  concentration. 
The  equation  for  the  friction  factor  that  follows  is  then 


1  _  (a 1 +A) 

/F "  v^ 


1-12 


For  a  given  flow  rate  and  polymer,  equation  I- 1 1  plots  as 


a  parallel  line  to  the  Newtonian  curve.  A  detailed  examination  of  the 


■ 

. 


-  8  - 


data  indicate  however  (especially  for  large  amounts  of  drag  reduction) 
that  the  parallel  lines  are  only  a  crude  representation  of  the  data. 
Alternatively  if  the  best  fit  straight  line  is  drawn  through  a  given 
velocity  profile  an  increase  of  the  slope  A  would  be  necessary. 

A  number  of  other  studies,  based  on  Pitot  tube  readings 
(49,  11,  P.6, 12)  have  subsequently  shown  velocity  profiles  which  are 
qualitatively  similar  to  Elata's. 

In  all  these,  owing  to  the  large  quantitative  uncertainty  in 
the  Pitot  readings,  it  is  impossible  to  determine  if  the  slope  is 
actually  different  or  if  modified  slope  simply  reflects  errors  in 
the  Pitot  readings.  It  has  been  shown  for  laminar  flows  (2  )  that  the 
Pitot  tube  measurements  in  a  viscoelastic  fluid  are  made  up  not  only 
of  the  kinetic  contribution  but  also  largely  influenced  by  normal 
stresses.  The  magnitude  of  the  normal  stress  contribution  depends 
on  the  local  shear  rate  or  radial  position  and  so  the  Pitot  readings 
may  not  even  determine  the  correct  shape  of  velocity  profile.  In 
turbulent  flows  the  flat  velocity  profile  in  the  core  region  would 
tend  to  hide  the  effect  of  shear  rate  but  would  be  serious  near  the 
wal  1 . 

Ni codemo  et  al  (29)  have  presented  velocity  profiles  using 
a  probe  calibrated  in  a  tow  tank.  It  was  decided  that  the  velocity 
profiles  could  not  be  interpreted  according  to  Bogue's  equation  1-5 
by  a  simple  increase  in  intercept.  Although  the  check  of  integrated 
and  measured  flow  rates  is  one  of  the  few  studies  showing  good  agree¬ 
ment,  this  check  is  a  very  insensitive  one  for  shape  of  profile. 

It  should  be  noted  their  calibration  does  not  indicate  the  magnitude 
of  any  contributions  owing  to  turbulence. 
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Similar  uncertainties  occur  in  the  case  of  heated  probles 
(25,13,23,33,40).  As  well  indicated  by  the  recent  study  of  Serth  and 
Kiser  (40)  errors  as  large  as  250%  in  determination  of  velocity  can 
occur.  Generally  these  studies  indicate  directly  there  is  an  undeter¬ 
mined  quantitative  uncertainty  in  all  reported  velocity  and  turbulence 
measurements  utilizing  either  pitot  or  heated  probes.  Qualitatively 
there  may  even  be  errors  in  shape  of  profile  owing  to  changes  in  probe 
response  with  radial  position. 

In  1969  Seyer  and  Metzner  (43)  presented  experimental 
velocity  profiles  obtained  from  a  photographic  method  which  it  is 
significant  to  note,  are  free  from  the  uncertainties  of  Pitot  or 
other  probes.  A  theoretical  analysis  similar  to  the  arguments  of 
Millikan  (27)  and  Dodge  (10)  was  used  to  show  that 

U+  =  A  In  Y+  +  B(e+) 

+  *?  1-13 

where  e  =  eU  Vv 

and  /j  =  A(l-02  In  (NRe/f)  +  (1  )2[B(e+)  -  A  In  2/2]  -  G  1-14 

where  ^  is  the  boundary  layer  thickness  (dimensionless)  and  G  is 
assumed  a  constant  equal  to  3.0.  Details  of  this  analysis  and  some 
modifications  will  be  presented  in  Cahpter  II. 

In  the  limit  for  C  and  going  to  zero  these  equations 
reduce  the  usual  equations  for  Newtonian  fluids  (1-1  and  1-6). 

For  drag  reducers  the  form  of  the  equations  is  the  same  as  Elata's 
with  the  exception  of  the  intercept  function. 


. 
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Seyer's  measurements  as  well  as  some  later  ones  in  a  square 
tube  by  Rudd  with  a  Laser  Doppler  Meter  (36),  are  in  excellent 
agreement  with  1-13.  In  each  of  these  studies  however  a  velocity 
profile  was  obtained  at  only  one  flow  rate  in  a  given  tube  for  turbu¬ 
lent  flow  conditions.  Therefore  they  do  not  provide  a  detailed 
verification  of  the  form  of  equations  1-13  and  1-14  nor  of  the  limits 
of  applicability. 

Seyer's  measurements  in  a  one-inch  tube  show  that  transitional 
flow  occurs  to  Reynolds  numbers  in  excess  of  10  for  what  have  been 
termed  concentrated  or  maximum  drag  reducers.  At  high  Reynolds  numbers 
the  function  factors  seem  to  approach  an  asymptote  defined  by  a  constant 
value  of  B(e+)  =  32.  Virk  et  al  (47)  have  published  a  number  of 
papers  which  consider  a  different  maximum  drag  reduction  asymptote. 

Based  on  data  from  a  number  of  sources  of  friction  measurements  it  was 
found  that 


1_ 

/f 


19  log  NRe  /F  -  32.4 


1-15 


Seyer's  friction  data  in  transitional  flow  follow  closely  to  the 
friction  factors  defined  by  equation  1-15.  The  corresponding  "ultimate 
velocity  profile"  is 


U+  =  11.7  In  Y+  -  17 


1-16 


Notice  this  is  a  much  higher  slope  than  for  Newtonian  fluids  and  has 
a  negative  intercept. 

The  transitional  flow  velocity  profiles  suggest  the 
interpretation  of  1-16  is  incorrect.  In  fact  in  the  final  limit  the 


' 
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profiles  would  be  laminar  over  the  entire  tube  cross  section.  Laminar 

+  + 

(parabolic)  profiles  on  U  -  Y  coordinates  are  easily  shown  to  be  a 
family  of  parallel  curves  depending  on  the  Reynolds  number.  For 
positions  near  the  center,  the  curves  can  be  approximated  by 
equations  like  1-16  with  large  slope  and  negative  intercept  as 
determined  by  the  Reynolds  number.  The  transition  profiles  on  these 
plots  are  qualitatively  similar.  Thus  Virk's  constants  may  simply 
reflect  the  transitional  nature  of  the  velocity  profiles.  In  the 
absence  of  defined  limits  of  transitional  flow  however,  little  more 
can  be  said  about  the  correct  interpretation  of  the  velocity  profiles 
for  transitional  flow.  Clearly  these  questions  are  of  prime  impor¬ 
tance  in  formulating  the  correct  similarity  laws  and  their  limits  of 
applicability. 
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1-3  PROPOSED  WORK 

In  this  work  it  is  proposed  to  extend  Seyer  and  Metzner's 
optical  method  of  measurements  in  order  to  obtain  instantaneous 
velocities,  both  axial  and  radial,  which  will  provide  a  detailed 
test  of  the  available  similarity  laws.  The  wall  region  as  well  as 
the  core  region  will  be  investigated  with  different  concentrations 
of  polymeric  solutions  for  different  pipe  sizes. 

The  statistical  properties  of  the  instantaneous  velocities 
will  be  used  to  comment  on  the  ’bursting  phenomena1,  and  any  change 
in  its  structure,  owing  to  the  presence  of  polymers. 
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THEORY 


II-l  PROBABLE  MECHANISM 

As  previously  mentioned,  the  recent  physical  interpretation 
of  turbulence,  presented  by  Corino  and  Brodkey  and  based  on  the 
observation  of  a  'bursting  phenomena1  present  in  the  wall  region  can 
serve  as  a  basis  to  analyse  the  data  collected  in  this  work.  Figure 
II-l  shows  schematically  a  lump  of  fluid  originating  in  the  main  flow 
and  entering  at  a  small  angle  into  the  boundary  layer  with  an  axial 
velocity  component  correspondi ng  to  the  average  velocity  of  its  origin. 
A  second  lump  of  fluid  located  in  the  boundary  layer  and  possessing  a 
lower  velocity  is  then  accelerated  by  the  excursion  and  shear  stresses 
start  to  build  until  one  or  more  ejections  occur.  Because  the  origin 
of  the  ejections  is  assumed  to  be  in  a  region  of  low  momentum 
situated  approximately  at  a  Y+  of  10,  the  axial  velocity  components 
of  the  ejections  are  expected  to  be  smaller  than  those  of  the  accel¬ 
erating  lumps.  Corino  and  Brodkey  observed  that  the  ejected  fluid, 
although  accelerated  in  the  axial  direction,  never  reached  the  axial 
component  of  the  main  flow  adjacent  to  the  boundary  layer. 


. 
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ILLUSTRATION  OF  BURSTING 
PROCESS 


wall 


FIGURE  II-1 
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Denoting  and  u^  as  the  velocity  components  of  the  lump 
of  fluid  entering  the  boundary  layer  and  u^  and  u  2  for  the  ejection 
velocity  components,  one  should  expect  to  find 

“xl  Cr)  >  “x2  (r)  "-1 

where 

7x  (r)  =  n  [=1(ux(r))i 

u  (r)  can  represent  the  time  average  velocity  at  a  fixed  radial  posi¬ 
tion  if  sufficient  readings  are  taken  (of  instantaneous  velocities) 
at  random  times.  Thus  a  histogram  of  the  axial  instantaneous 
velocities,  at  a  fixed  radial  position,  should  be  composed  of  two 
distributions.  Using  a  flash  photolysis  method,  Popovich  and  Hummel 
(33  ,34)  measured  velocity  profiles  which  showed  a  binodal  distribution 
of  the  axial  instantaneous  velocities  in  the  boundary  layer  with  the 
peak  of  lower  axial  velocities  becoming  more  pronounced  at  positions 
nearer  the  wall.  This  simply  shows  the  influence  of  the  velocities 
of  fluid  originating  near  the  wall  dominates  the  flow  in  that  region. 

The  ejection  phenomena  approach  taken  in  this  work  can 
explain,  in  more  detail,  the  presence  of  dual  distributions  than  the 
Popovich  analysis,  which  is  based  on  elements  of  fluid  being  directed 
toward  the  wall  and  others  directed  toward  the  center. 

Each  positive  (directed  toward  the  wall)  radial  velocity  is 
associated  with  a  u  |  while  a  negative  radial  velocity  is  associated 
with  ux2.  The  velocity  histograms  can  then  be  split  according  to 
whether  the  instantaneous  radial  velocity  is  positive  or  negative. 
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According  to  the  bursting  analysis,  instantaneous  negative  radial 
velocities  should  be  larger  than  the  positive  components.  Thus  to 
obtain  a  zero  time  averaged  radial  velocity  (the  sum  of  all  observed 
components)  the  positive  components  must  be  observed  more  often  than 
the  larger  negative  components.  Defining 

RQ  =  n+/n"  1 1-2 

where 

n+  is  the  observed  number  of  u^ 

n  is  the  observed  number  of  u  0 

r2 

then  Rq  is  expected  to  be  greater  than  unity  in  the  boundary  layer  and 
approximately  equal  to  unity  in  the  core  region.  Also  if  the  velocity 
of  ejection  is  dependent  on  radial  position,  as  Corino  and  Brodkey 
observed,  an  acceleration  of  the  fluid  during  the  trajectory  from  its 
origin  to  the  radial  position  where  it  is  broken  up  by  the  main  flow, 
implies  the  presence  of  a  peak  in  the  plot  of  Rq  versus  dimensionless 
radial  position.  This  peak  should  occur  in  the  vicinity  of  the  edge 
of  the  boundary  layer. 

In  summary,  an  analysis  of  the  histograms  of  axial  and  radial 
components  of  the  instantaneous  velocities  at  a  radial  position  in  the 
boundary  layer  should  provide  information  about  the  mechanism  of 
ejection.  Also  it  should  explain  the  presence  of  the  binodal  distri¬ 
butions  in  the  boundary  layer,  if  these  distributions  are  split 
according  to  the  direction  of  the  radial  component.  In  the  case  of 
drag  reducers  the  available  data,  although  in  doubt  quantitatively, 
show  that  the  turbulence  is  similar  to  Newtonian  turbulence.  The 
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drag  reduction,  therefore,  is  probably  a  result  of  modifications  to  the 
ejection  process.  For  example,  if  the  boundary  layer  is  thickened 
appreciably,  then  the  peak  in  Rq  ought  to  appear  at  greater  values  of 
Y+  than  in  the  Newtonian  case. 

In  the  ejection  process,  earlier  studies  (37,43)  suggest  the 
increased  resistance  to  stretching  causes  the  radial  fluctuations  to 
decrease.  Although  in  the  work  of  Seyer  and  Metzner,  where  the  ideal 
case  of  steady  stretching  was  treated,  the  same  conclusion,  that  is  an 
increased  resistance  to  stretching,  follows  from  a  consideration  of 
the  transient  motion  (8).  In  Figure  1 1 - 1  the  resistance  to  sudden 
deformations  would 

(i)  decrease  u  at  a  given  point, 

r2 

(ii)  not  affect  appreciably  u  and  u  ,  because  the  direction  of 

rl  X1 

these  larger  lumps  of  fluid  are  nearly  parallel  to  the  flow 

direction  (no  stretching  in  the  axial  direction), 

(iii)  implies  a  lower  radial  variance  from  a  combination  of  (i )  and 

(ii)  (u  not  affected,  and  u  decreased). 
rl  r2 

One  cannot  predict  the  change  in  u  but  if  a  decrease  results,  as 

x2 

might  be  expected  if  the  sublayer  thickens  appreciably,  then  higher 
axial  intensity  would  result. 

Similarly,  however,  the  data  of  Meak  and  Baer  (24)  and 
Fortuna  and  Hanratty  (15),  show  the  frequency  of  occurrence  of  the  ejec¬ 
tion  decreases  in  polymers.  This  would  also  decrease  the  radial  intensity. 
Unfortunately,  the  method  used  in  this  work  cannot  distinguish  directly 
which  cause  is  effective  in  altering  the  intensity. 
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1 1-2  SIMILARITY  LAWS 
A)  Velocity  Profiles 

Assuming  the  relaxation  time  as  the  characteristic  time  scale 
of  viscoelastic  solutions,  Seyer  and  Metzner  (43)  introduced  the 
relaxation  time  as  an  additional  parameter  in  the  wall  region  where 
the  stretching  is  assumed  to  occur.  They  performed  an  analysis, 
similar  to  Millikan  (27),  to  obtain  the  logarithmic  velocity  profile. 
The  velocity  in  the  boundary  layer  can  be  expressed  as 

ux  -  f2  (v  y>  p.  n>  s)  n'3 

and  in  the  core  region  as 

<Vm  -  “x  =  *3  (V  R>  y*  p>  n"4 

so  that  for  the  entire  section  of  the  pipe 

ux  =  F  (tw,  R,  y,  p,  u,  0)  I 1-5 

Notice  that  the  fluid  parameters  should  be  estimated  at  the  wall  shear 
stress  value.  Introducing  the  following  dimensionless  group 

Z  =  RU*/v  1 1 -6 

and  using  those  already  defined  (in  Chapter  I),  the  following  equations 
for  the  core  and  the  wall  regions  can  be  found  from  a  dimensional 
analysis 


-  u. 


<u  > 
x 


-  A  In  £  +  C3 


m 


1 1 -7 


. 
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and 

uv  /U*  =  A  In  Z  +  B(e+)  1 1-8 

x  m 

when  combined,  the  expression  for  the  overlapping  region  is  obtained 

U+  =  A  In  Y+  +  B(e+)  -  C3  1 1-9 

Extensive  study  was  performed  by  Bogue  (5)  to  find  an  expres¬ 
sion  for  the  factor.  must  be  a  function  of  £  to  permit  equation 
1 1-7  to  hold  for  the  entire  region.  It  must  be  zero  at  the  centerline 
and  contribute  approximately  10%  to  values  of  the  velocity  defect  for 
£  smaller  than  0.3  as  shown  by  Bogue.  In  equation  I I -9  the  correction 
function  may  be  written  in  the  form 

_  .  (S-0.8)2 

C3(£,f)  =  0.05/y  e  0,‘5  11-10 

so  that  the  corrected  velocity  profile  becomes 

U+  =  A  In  Y+  +  B(e+)  +  C3(?,f)  11-11 

and  where  for  Newtonian  fluids  A  =  2.46  and  B(0)  =  5.6.  Assuming, 
as  an  approximation,  a  linear  profile  in  the  boundary  layer  and  neg¬ 
lecting  C3(£,f)  for  small  £,  Equation  11-11  becomes 

Yj  =  A  In  Y|  +  B(e+) 


11-12 
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and  serves  to  define  the  boundary  layer  thickness. 

Figure  1 1 -2  shows  the  Cc  »f )  function  increasing  with 
increasing  wall  distance,  £.  It  can  be  observed  that  the  function 
increases  in  value  as  the  friction  factor  decreases,  indicating  that, 
if  applicable,  the  correction  becomes  more  important  for  drag 
reducers.  Because  of  the  passive  role  of  the  core  region  when  drag 
reduction  occurs,  as  shown  experimentally  by  Wells  and  Spangler  (48) 
and  suggested  by  the  bursting  arguments,  the  same  function  may  apply 
for  viscoelastic  fluids. 

B)  Friction  Factor 

In  order  to  express  the  friction  factor  defined  in  Equation 
1-4  ,  the  bulk  velocity  is  obtained  from  the  integrated  profile  which 
can  be  divided  into  the  boundary  layer  region  and  the  core  region 


<ux> 


R 


r  u  (r)  dr  + 


R 


r  uxc(r)dr 
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In  terms  of  the  dimensionless  radial  position 


<u  >  =  2 

x 


1  (l-?)uxLU)d?  +  2 


Ci-5)uxcCe)de 


11-14 


where 


51  =  6/R 

Assuming  a  linear  profile  in  the  wall  region  and  using  the 
defect  law  for  the  core  region,  one  can  obtain 
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FIGURE  II - 2 


BOGUE'S  CORRECTION  FUNCTION 
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22  - 


<V 


RU 

v 


‘  2  3" 

ii  ii 
2  '  3 


<u  > 

X 


-  (1-c,) 


o  * 
^  -  U 


rl 


(l“5)f3U)d5  11-15 


1 


Combining  Equations  1-4  and  11-15  and  grouping  terms,  the  following 
expression  is  obtained 

fj  =  A(l-5l)2  ln(NRe/f)  +  (l-C,)2[B(e+)  -  A  In  2/2] 


G(Cr  ^  ,  f) 


11-16 


which  is  the  same  as  Seyer's  equation  (41  )except  that  G  is  now  also  a 
function  of  f.  The  function  G  can  be  expressed  as 

■1  fSi 


G  =  2  I  (l-5)f,U)dC  -  2  1  (l-Of3U)de  -  52  ^ 

J  0  0 


11-17 


where 
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(l-?)fo(c)dc  =  -2A 

In  c  dc  - 

C  In  C  d c 
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C3(c,f)d?  -  2 


c3U,fk  dc 


In  Equation  11-18  with  A  =  2.46 


2A 
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The  curves  of  Figure  1 1-2  are  used  to  graphically  integrate  the  remain¬ 
ing  terms  of  Equation  11-18  involving  C3(c,f)  giving  a  value  of  0.40 
for  the  case  of  a  very  low  friction  factor  of  0.003  where  C3U,f)  is 
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largest.  Because  the  correction  factor,  ,f)  is  zero  for  small 
values  of  the  remaining  terms  of  Equation  11-17  are  not  influenced 
by  the  correction  factor.  Therefore,  since  G  is  a  minor  term  in 
Equation  11-16,  the  dependence  of  G  on  f,  as  determined  by  Cg(j;,f), 
is  negligible. 

Figure  1 1-3  shows  the  B  function  evaluated  from  friction 
factors  and  velocity  profiles  by  Seyer  (41)  for  solutions  of  Separan 
AP-30  for  1"  and  2"  pipes  at  three  different  concentrations.  From 
Figure  1 1 - 3  it  can  be  observed  that  the  function  reaches  an  asymptotic 

-f 

value  of  approximately  32.2  for  e  greater  than  20.  If  this  function 

is  unique  for  all  viscoelastic  fluids  and  pipe  sizes,  then  maximum 

+ 

drag  reduction  will  be  reached  whenever  e  is  greater  than  20.  The 
only  limitation  on  the  use  of  Equation  1 1 - 1 1  is  the  dividing  point 
between  turbulent  and  transitional  flow. 

The  limits  of  transitional  flow  must  be  related  to  the 
thickening  of  the  sublayer.  The  definition  of  sublayer  thickness 
given  by  Equation  11-12  shows  how  it  depends  directly  on  the  elas¬ 
ticity  of  the  system.  Clearly  the  sublayer  must  be  less  than  the  tube 
radius  for  turbulent  flow  to  occur,  but  the  question  of  an  exact 
value  needs  further  investigation.  In  the  case  of  Newtonian  fluids 
the  maximum  sublayer  thickness  is  approximately  £  -  0.1  (at  N^e  of  2100). 


1 


FIGURE  II- 3 


-  24  - 


o  o  o  O  o  o 

LO  XT  CO  CNJ  t— 

m 


-  25  - 


EXPERIMENTAL 

An  optical  method  of  measurement  was  used  to  obtain  streak 
photographs  of  small  air  bubbles  in  the  water  and  polymeric  solutions 
in  order  to  measure  instantaneous  velocities  in  the  boundary  layer 
and  core  regions.  As  mentioned  before,  this  method  does  not  suffer 
from  the  uncertainties  inherent  in  pitot  tube  and  heated  probe 
measurements . 

III-l  FLOW  SYSTEM 

The  flow  system,  as  shown  schematically  on  Figure  III-l  ,  con¬ 
sisted  of  a  closed  circuit  in  which  the  fluid  could  be  pumped  through 
a  test  section  and  returned  to  a  300  gal.  stainless  steel  storage 
tank.  A  variable  speed  Lightnin  mixer  and  a  coil  were  attached  to 
the  tank  to  regulate  the  level  of  degradation  and  the  operational 
temperature  of  the  Separan  AP-30  solutions.  Two  Moyno  positive  dis¬ 
placement  pumps,  each  with  a  maximum  delivery  capacity  of  200  gal /min 
permitted  measurement  of  pressure  drop  over  an  extended  range  of  flow 
rates.  A  2"  I.D.  Foxboro  magnetic  flowmeter  recorded  the  flow  rate 
passing  through  the  test  section. 

A  variable  pressure  Desurger  following  the  pump  exit  damped 
the  flow  fluctuations  caused  by  the  positive  action  of  these  pumps. 
Nominally  the  desurger  pressure  was  set  for  the  photographic  runs  at 
75%  of  the  pump  discharge  pressure.  A  calming  section  of  6"  I.D.  and 
6'  long  containing  a  baffle  plate  was  inserted  before  the  test  section 
in  order  to  fulfill  four  requirements: 

i)  injection  point  of  the  tiny  air  bubbles  used  as  tracers. 
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FIGURE  III-1 
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ii)  trapping  of  the  dirt  particles  in  the  system, 

iii)  trapping  of  large  air  bubbles, 

iv)  calmed  the  fluid  and  uniformly  distributed  the  tracers  in  the 
fl  uid. 

Air  bubbles  injected  into  the  fluid  flow  were  used  as 
tracers.  The  injection  point  was  situated  in  the  calming  section  in 
front  of  the  baffle  plate  (see  Figure  III -2 ) .  The  number  and  size 
of  the  air  bubbles  were  adjusted  by  a  needle  valve  and  a  porous 
chamber  located  inside  the  calming  section.  This  variable-position 
porous  chamber  allowed  the  bubbles  to  be  mixed  completely  in  the  flow 
by  locating  it  in  the  large  eddies  formed  at  the  front  of  the  plate. 
Additionally,  since  the  velocity  of  the  flow  was  very  small  in  this 
6"  section,  large  air  bubbles  migrated  to  the  top  where  they  could 
be  purged  from  the  system.  Simultaneously,  dirt  and  rust  particles, 
which  caused  a  loss  of  contrast  in  photographs,  could  be  collected 
at  the  bottom. 

These  variables  were  adjusted  to  obtain  an  average  diameter 
of  the  bubbles  of  approximately  0.002".  The  maximum  observed  diameter 
of  photographed  bubbles  was  0.0041"  in  diameter,  which  is  smaller  than 
the  smallest  turbulence  scale.  In  all  cases,  however,  readings  of  the 
large  bubbles  (0.004")  were  rejected. 

The  maximum  percentage  volume  of  air  in  the  fluid  was  cal¬ 
culated  using  photographs  obtained  at  an  exposure  time  of  1/1000 
second  for  run  14.  On  the  average,  18  bubbles  appeared  per  frame, 
between  the  wall  and  the  center  of  the  pipe.  A  percentage  of  air 
by  volume,  the  maximum  relative  to  the  size  of  the  bubbles  and  to  the 
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estimated  volume  of  fluid,  was  then  calculated,  using  0.004"  as  the 
diameter  of  the  bubbles.  Thus  all  bubbles  are  assumed  to  be  of 
maximum  size  in  the  photographed  volume  of 

Vfluid  =  °“5"  x  °-752"  x  0-0592"  =  0.0222  in3 
The  volume  of  air  is 

V  .  =  18(0.524  x  0.0043)  =  604  x  10"9  in3 

a  i  r 

and  therefore 

%  volume  of  air  =  (604  x  10"9)  x  100/0.0222  =  0.00272 

which  is  considered  negligible.  It  is  interesting  to  note  that, 
at  the  position  where  photographs  were  obtained,  the  air  bubbles 
were  found  to  be  unevenly  distributed  in  the  illuminated  cross 
section  of  the  pipe.  A  maximum  concentration  occurred  at  a  radial 
position  r/R  =  0.890  (see  Appendix  G  for  the  experimental  distribu¬ 
tion  curves  of  radial  concentration). 

Two  different  systems  of  test  section  were  designed  in 
order  to  study  the  diameter  effect  present  with  drag  reduction.  The 
first  system  consisted  of  precision  bore  pyrex  tubes  of  1"  inside 
diameter  and  the  second  system  consisted  of  standard  plexiglass 
pipes  of  2.75"  inside  diameter  (dimensions  of  these  systems  can  be 
found  in  Table  1 1 1 -1 ) .  It  may  be  noted  in  Table  1 1 1 - 1  that,  in  the 
case  of  the  second  system,  part  of  the  1"  tube  was  used  and  connected 
as  a  bypass,  in  order  to  measure  pressure  drops  and  determine  the 
level  of  degradation  of  the  solutions  being  used  in  the  2.75"  pipe. 
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TABLE  III-l 

TEST  SECTIONS 


System  I  System  II 


Nominal  diameter  (inch) 

1 

1 

2.75 

Materi al 

pyrex  bore 

pyrex  bore 

plexiglass 

Inside  diameter  (inch) 

1 .0  +  0,001 5 

1 .0  +  0,0015 

2.75  +  0.005 

Total  length  (ft) 

28 

21 

36 

Entry  length  to 
optical  box  L/D 

187 

-- 

114 

Entry  length  to  first 
pressure  taps  L/D 

54 

66 

60 

Pressure  taps  -hole 
diameter  (inch) 

0.025  +  0.005 

0.025  +  0.005 

0.020  +  0.005 

Length  between 
pressure  taps  (in  ft) 

OM4) 

22 

23.60 

(2 )- (4) 

16 

12 

17.50 

(3 )- (4) 

11 

6 

11.69 

• 

oa.ts 
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Four  sections  of  precision  bore  1“  tubes  7  feet  in  length 
were  connected  and  clamped  on  a  lathe  bed  in  order  to  cut  down  vibra¬ 
tions  originating  mainly  from  the  pump.  Great  care  was  taken  to 
ensure  that  the  upstream  joints  presented  no  discontinuity  to  disturb 
the  flow.  This  was  done  by  inserting  a  Teflon  plug,  0.995"  in 
diameter  and  3"  in  length,  inside  the  tubes  at  the  joint  before  con¬ 
necting  two  sections  together.  Four  pressure  taps  were  used  and 
connected  to  two  60"  manometers  containing  Mercury  and  0-Dichloro 
Benzene  of  density  1.305. 

Similarly,  in  replacing  the  1"  system,  the  2.75"  sections 
were  carefully  jointed  and  clamped  to  the  lathe  bed.  Three  pressure 
taps  were  connected  to  the  same  manometers.  The  justification  for 
the  many  pressure  drop  measurements  along  the  pipe  was  to  check  the 
reproducibility  of  the  measurements  and  also  to  verify  that  the  large 
entry  length  (  6)  had  been  exceeded  at  the  optical  section. 


.£  - 
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1 1 1-2  EXPERIMENTAL  FLUIDS 

The  drag  reducing  fluid  chosen  for  the  entire  work  was  a 
polyacrylamide  product  of  Dow  Chemical  Co.,  named  Separan  AP-30. 

This  granular  white  powder,  of  an  estimated  molecular  weight  between 
2  and  3  million,  possesses  elastic  properties  when  in  solution  with 
water.  It  was  preferred  to  other  polymers  because  of  its  known  level 
of  drag  reduction  which  can  be  found  in  previous  works  (30  ,43  ). 

Two  different  weight  concentrations  were  used:  the  0.01% 

(dilute)  and  the  0.10%  (concentrated).  Physical  properties  of  these 

solutions  are  presented  in  Appendix  A.  The  viscous  properties  of  the 

0.01%  solutions  showed  that  their  behavior  can  be  assumed  to  be 

-5  2 

Newtonian  with  a  kinematic  viscosity  of  1.2  x  10  ft  /sec  at  room 
temperature. 

To  ensure  the  same  degradation  level  for  all  of  the  final 
data,  a  new  solution  was  prepared  for  each  run  and  the  following  pro¬ 
cedure  was  used.  A  mechanical  degradation  by  shearing  in  the  1" 
pipe  at  a  high  flow  rate  was  followed  until  the  recorded  pressure 
drop  measurements  reached  the  drag  reduction  level  desired.  The 
desired  level  was  determined  arbitrarily  by  a  master  curve  of  pres¬ 
sure  drop  versus  flow  rate  in  the  1"  pipe  which  had  been  determined 
so  that  enough  drag  reduction  could  be  obtained  in  the  2.75"  pipe. 
Then  the  photographic  run  was  started  with  a  check  of  pressure  drop 
before,  in  the  middle,  and  at  the  end  to  assure  that  negligible 
degradation  occurred.  A  slightly  different  procedure  was  followed 
for  the  determination  of  the  level  of  degradation  of  the  solutions 
used  in  the  2.75"  pipe.  Here  the  solutions  were  degraded  in  the  l1 
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pipe  to  the  desired  level,  the  flow  was  then  diverted  to  the  2.75"  pipe 
for  readings  of  pressure  drop  and  photographs  taken.  This  ensured  that 
drag  reduction  and  velocity  measurements  in  the  2.75"  pipe  would  be 
at  the  same  level  of  degradation  as  measurements  in  the  1"  pipe. 
Friction  measurements  were  also  checked  at  the  end  of  a  photographic 
run  in  the  2.75"  pipe  by  additional  measurements  in  the  1"  pipe. 
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1 1 1-3  OPTICAL  SYSTEM 

The  main  features  of  the  optical  system  shown  in  Figure  1 1 1-3 
are 

i)  a  light  source  which  illuminated  a  thin  horizontal  portion  of 
the  pipe  cross  section, 

ii)  a  chopping  disk  for  the  determination  of  the  time  scale, 

i i i )  an  optical  box , 

iv)  a  camera  assembly. 

The  instantaneous  velocities  of  the  fluid  in  the  pipe  were 
determined  by  photographing  the  air  bubbles  sweeping  through  the  narrow 
light  beam  located  at  the  center  of  the  pipe.  A  three  phase  'Christie' 
transformer  of  D.C.  maximum  output  of  3300  watts  at  100  amps  operated 
a  'Christie'  Xenon  arc-lamp.  The  converging  mirrors  in  the  lamp¬ 
housing,  No.  BSF  50,  were  positioned  to  direct  the  light  beam  through 
the  view  section  of  the  pipe.  An  adjustable  slit  was  placed  immediately 
in  front  of  the  lamphousing  lens  to  produce  a  narrow  beam,  and  converg¬ 
ing  lenses  used  after  to  adjust  the  light  beam  to  a  thickness  smaller 
than  the  slits  on  the  chopping  disk. 

Aluminum  disks  of  15"  in  diameter  and  1/8"  thick  were 
assembled  on  three  phase  motors  to  obtain  constant  rotational  speed. 

In  order  to  interrupt  the  light  and  provide  streaks  of  known  duration, 
each  wheel  was  divided  into  an  equal  number  of  radial  slits  and  spokes. 
The  rotational  speed  of  the  motors  and  the  number  of  slits  on  the 
wheels  were  determined  in  order  to  obtain  time  scales  smaller  than  the 
reciprocal  of  the  estimated  frequencies  of  turbulence.  Table  III -2 
represents  the  different  combinations  used  for  the  various  runs. 
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FIGURE  III-3 


OPTICAL  SYSTEM 
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These  time  scales  also  provided  convenient  streak  lengths  on  the  films. 


TABLE  1 1 1-2 


Motor 

rpm 

Number  of  Divisions 
on  the  wheel 

Time  Scale 
sec. 

Run  Number 

1790 

40 

0.001675 

4, 5, 6, 7 

3577 

40 

0.000837 

9,10 

3575 

60 

0.00056 

11  ,12,13,14 

1790 

60 

0.001118 

15,16,17 

The  chopping  disks  were  placed  perpendi cul ar  to  the  light 
beam  in  front  of  the  view  section  with  the  axis  of  the  motor  in  the 
plane  of  the  light  slit.  Optical  boxes  placed  around  the  test  sections 
were  slightly  different  in  construction.  For  the  1"  pipe,  the  view 
section  consisted  of  a  plexiglass  box  filled  with  a  mineral  oil  of 
approximately  the  same  index  of  refraction  as  glass.  The  second  type 
consisted  of  a  polished  cylinder  of  2.75"  I.D.  in  a  cubic  piece  of 
plexiglass  which  was  inserted  between  two  sections  of  pipe.  The 
upstream  joint  was  then  polished  to  remove  any  discontinuity. 

Figure  1 1 1-4  shows  a  scheme  of  the  optical  box  used  for  the 
2.75"  pipe.  A  slit  of  0.0592"  in  height  and  2"  in  width  was  placed  on 
the  view  section  to  assure  the  illumination  of  a  volume  of  fluid  of 
nearly  the  same  X-section.  The  outside  of  the  optical  box  was  painted 
black  in  order  to  cut  down  possible  reflections. 

A  needle  of  0.020"  in  diameter,  inserted  in  the  test  section 
through  a  partially  sealed  hole,  could  be  moved  across  the  illuminated 
portion  of  the  fluid.  A  cylindrical  guide  was  attached  to  the  rear 
portion  of  the  view  section  to  properly  align  the  needle.  This 
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needle  served  to  obtain  the  radial  calibration  curves  (see  Appendix  B) 

A  Miranda  35  mm  camera,  located  perpendicular  to  the  light 
beam,  was  attached  to  Setco  Building  Blocks  equipped  with  micrometer 
handwheels  graduated  in  0.001"  increments.  Real  magnifications  of 
approximately  8X  and  2X  were  obtained  from  different  combination  of 
extension  tubes  and  lenses.  A  1:1.9/50  mm  lens  with  a  depth  of  field 
of  0.077"  was  used  for  the  runs  in  the  1"  pipe  while  a  1:3.5/80  mm 
lens  with  a  depth  of  field  of  0.65"  was  used  in  the  2.75"  pipe. 
Thirty-five  mm  Tri-X  films  with  ASA  speed  of  400,  were  used  in  this 
study. 

The  best  exposure  times  and  shutter  openings  were  deter¬ 
mined  by  a  trial  and  error  procedure; 

a)  obtaining  a  maximum  number  of  streaks  per  frame  (one  streak  per 
radial  position  if  possible)  to  cut  down  the  labor  of  analysing  the 
films  and  to  economize  on  film.  This  necessitated  a  long  exposure 
time. 

b)  alternatively,  a  long  exposure  time  meant  that  the  illuminated 
curvature  of  the  wall  would  appear  very  dark  and  would  mask  the 
streaks  near  the  wall.  This  necessitated  closing  the  shutter 
opening . 

c)  finally  the  air  bubbles  used  as  tracers  had  to  appear  very  sharp 
on  the  films. 

Experimentally,  it  was  found  that  an  exposure  time  of  1/15 
of  a  second  with  a  shutter  opening  of  f:4  gave  the  best  results. 
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III -4  OPERATIONAL  PROCEDURE 

After  a  solution  was  found  to  be  at  the  desired  level  of 
degradation,  it  was  pumped  through  the  system  at  a  very  high  flow  rate 
for  a  few  minutes  to  enable  the  air  to  escape  from  the  apparatus.  The 
light  source  was  turned  on  and  the  optical  system  adjusted.  Then  the 
flow  rate  was  adjusted  to  the  desired  value  and  the  tracers  were 
introduced  into  the  system.  When  the  flow  rate  and  the  other  variables 
of  the  run  were  found  to  be  steady  (as  indicated  by  the  magnetic  flow 
meter,  the  thermometer,  etc.)  photographs  were  taken.  The  photographic 
runs  took  3  to  4  hours  for  the  first  few  runs,  but  only  between  1  and 
1.5  hours  for  all  other  runs  due  to  the  amelioration  of  the  procedure 
to  obtain  calibration  curves  (see  section  on  Analysis  of  Streak 
Photographs) . 


' 
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1 1 1-5  STREAK  PHOTOGRAPHS 

Typical  streak  photographs  are  shown  on  Figures  1 1 1-5  and 
1 1 1 -6  for  the  case  of  four  different  runs  at  four  different  magnifi¬ 
cations.  Figure  1 1 1 - 5a  indicates  an  instant  in  time  in  which  all  of 
the  streaks  are  parallel  to  the  wall.  This  streak  pattern  is  observed 
on  a  portion  of  the  photographs  during  the  photographic  runs  of  the 
concentrated  solutions  and  indicates  a  patch  of  laminar  flow.  Other 
photographs  which  show  large  radial  components  indicate  a  patch  of 
turbulent  fluid.  Identical  streak  photographs  can  be  found  in  Seyer's 
work  (43).  In  contrast,  Figure  1 1 1 - 5b  shows  the  presence  of  radial 
components  and  represents  a  fully  turbulent  flow.  Those  radial 
components  as  it  is  emphasized  from  an  enlarged  view,  as  shown  on 
Figure  III -6a  for  a  water  (run  II),  were  observed  on  every  photograph. 

Figure  1 1 1 -6b  shows  three  streak  patterns  located  in  an 
illuminated  region  (from  illumination  of  the  curved  inside  wall  of 
0.059l:  width),  to  indicate  how  close  to  the  wall  photographs  of 
tracers  could  be  obtained  at  high  magnification.  The  closest  streak 
pattern  is  located  at  radial  position  approximately  200  microns  from 


the  wal 1 . 
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FIGURE  5 


b  RUN  16  ( Ux)  -  5  09  FJ7SEC 


FIGURE  6 


a  RUN  II  (  0^=13.67  Ft/SEC 


b  RUN  t3/Ux)  =  21.62  FT/SEC 


-  43  - 


1 1 1-6  ANALYSIS  OF  STREAK  PHOTOGRAPHS 

All  photographs  were  analysed  using  a  'digitizer'  which 
determines  the  X-Y  coordinates  of  a  picture  projected  onto  a  glass 
screen.  The  projection  provides  an  additional  magnification  of 
approximately  20X.  The  'digitizer'  automatically  provided  an  IBM 
punched  card  containing  the  X-Y  coordinates,  in  volts,  of  points 
relative  to  the  'digitizer'  frame.  The  X-Y  frame  of  the  digitizer 
(22.9"  x  21.7"  in  full  scale)  consisted  of  a  cross  hair  device  con¬ 
nected  to  two  variable  resistances.  A  power  supply  coupled  to  the 
variable  resistances  provided  a  reading  of  10.000  ±  0.005  volts  for 
the  full  trajectory  of  the  axis.  Six  voltages  were  then  punched  on 
a  card  and  two  extra  numbers  were  added  in  order  to  define  the  number 
of  streaks-spokes  read  and  the  radial  position. 

Figure  III -7  shows,  schematically,  a  projected  photograph 
containing  a  single  streak  pattern.  The  large  arrows  represent  the 
interrupted  streak  and  its  direction  for  a  period  corresponding  to 
the  time  taken  by  three  spokes  of  the  timing  wheel  to  cut  the  light. 

The  small  lines  represent  the  axial  and  radial  components  of  the 
streak  and  spoke  for  one  time  period. 

In  order  to  determine  the  streak  length  in  the  axial  and 
radial  directions  for  a  known  time  period,  readings  of  four  coordinates 
relative  to  the  fixed  frame  of  the  digitizer  were  needed.  Points  1 
and  2  represent  the  beginning  and  end  of  a  streak-spoke,  respectively, 
and  points  3  and  4  locate  the  line  determining  the  pipe  wall  and  con¬ 
sequently  the  axial  direction  in  the  pipe.  With  these  readings  and 
appropriate  calibrations  which  convert  the  voltages  to  actual  distances, 
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FIGURE  III -7 


wall 
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the  resolution  of  a  streak  into  a  radial  and  an  axial  component  follows 
di rectly . 

The  dashed  lines  represent  the  band  defining  a  radial  posi¬ 
tion  where  the  streaks  were  classified  as  acceptable.  The  true  radial 
position  for  each  acceptable  streak  was  taken  to  be  the  location  of  the 
center  of  the  band.  In  the  core  region,  where  the  velocity  and  inten¬ 
sity  gradients  are  small,  the  width  of  the  band  is  of  little  consequence. 
However,  this  is  not  true  very  near  the  wall.  For  the  average  velocity 
near  the  wall,  since  the  profile  is  essentially  linear  over  the  width 
of  the  band,  the  average  is  not  changedo  However,  the  calculated 
intensity  will  be  slightly  higher  than  the  actual  intensity. 

Generally,  if  two  or  more  correlated  streak  patterns  were 
observed  in  the  band,  only  one  reading  of  velocity  was  taken.  This 
tends  to  eliminate  the  bias  of  too  many  readings  of  the  same  instant¬ 
aneous  velocity. 

The  reproducibility  of  readings  and  mechanical  errors  involved 
in  the  use  of  the  'digitizer'  are  presented  in  Appendix  C.  The  instant¬ 
aneous  velocities  were  calculated  on  the  IBM  1800  facilities  in  this 
department.  Programs  used  in  this  work  are  presented  in  Appendix  J 
which  can  be  found  in  Book  II. 

To  determine  the  axial  magnification,  pictures  of  a  Keuffel 
and  Esser  ruler  were  taken  outside  the  tube  at  identical  camera  focus 
(with  respect  to  the  photographic  runs)  and  projected  on  the  'digitizer' 
screen  where  measurements  of  the  scale  were  obtained.  Trial  runs  have 
shown  that  identical  magnifications  were  obtained  with  the  ruler 
inside  and  outside  of  the  pipe  (see  Table  B - II). 
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The  technique  developed  to  find  the  radial  magnification  and 
the  radial  positions  was  complicated  by  the  non-linearity  of  the  curves 
due  to  wall  curvature.  A  needle  of  0.020"  in  diameter  was  positioned 
perpendicular  to  the  camera  and  located  at  the  center  of  the  pipe  or 
light  beam  (see  Figure  1 1 1-4).  This  needle  was  attached  to  a  micro¬ 
meter  head  with  a  smallest  division  of  2  microns  and  mounted  on  a 
sliding  block.  The  needle  was  then  moved  to  the  wall  of  the  pipe 
(where  it  could  be  photographed)  and  a  sequence  of  pictures  were  taken 
as  the  needle  was  retrieved,  one  position  at  a  time.  The  fine  adjust¬ 
ment  micrometer  head  provides  accurate  positions  in  the  wall  region 
and  the  handwheels,  graduate  in  0.001"  increments,  provided  subsequent 
measurements  to  the  center  of  the  pipe.  For  most  of  the  runs,  this  set 
of  pictures  was  checked  by  a  second  set  which  was  obtained  with  the 
needle  moving  from  the  center  of  the  pipe  toward  the  wall.  Appendix  B 
presents  typical  data  for  the  calibration  curves  obtained  and  shows  a 
typical  example  of  the  third  order  polynomial  radial  calibration  curves. 
The  knowledge  of  both  axial  and  radial  magnifications  and  the  radial 
positions  permitted  the  voltages  to  be  transformed  into  instantaneous 
axial  and  radial  distances,  whereby  ax  and  Ar,  could  be  calculated. 

Knowing  the  time  scale  and  length  scale  for  each  streak- 
spoke,  the  axial  and  radial  instantaneous  velocities  could  be  found 

ux(r)  =  ax/T 

ur(r)  =  Ar/T 
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where  T  =  60/ (M )  (W)  —  (sec) 


M  =  number  of  slit-spoke  on  the  chopping  wheel 

All  the  instantaneous  velocities  are  presented  in  Appendix  K,  which 
can  be  found  in  Book  II. 

The  mean  velocities  could  then  be  calculated,  at  a  radial 
position,  from 
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x  n 


and  the  Root  Mean  Square  values  calculated  from 
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RESULTS  AND  DISCUSSION 

IV-1  PRESSURE  DROP  MEASUREMENTS 

Figure  IV-1  represents  the  friction  factor  -  generalized 
Reynolds  number  data  obtained  for  both  water  and  polymeric  solutions 
in  the  1"  and  2.75"  diameter  pipes.  These  data  are  tabulated  in 
Appendix  I. 

The  experimental  values  of  the  friction  factor  for  water 
were  obtained  in  order  to  check  the  experimental  equipment  as  well  as 
the  Newtonian  form  of  the  friction  similarity  law.  All  the  Newtonian 
data  representing  the  curve  labelled  water,  in  Figure  IV-1,  are 
obviously  in  close  agreement  with  the  accepted  Newtonian  curve. 

In  Figure  IV-1  considering  the  0.01%  solution  first,  the 
data  represent  a  certain  level  of  degradation  of  the  polymer.  As 
shown  by  comparison  of  friction  factors,  the  amount  of  degradation  in 
this  work  is  somewhat  less  than  the  amount  of  degradation  in  Seyer 
and  Metzner's  work  with  equivalent  solutions.  The  somewhat  'fresher' 
solutions  used  were  necessary  in  order  to  get  significant  drag 
reduction  in  the  2.75"  tube.  If  the  solutions  were  too  fresh  however, 
it  was  found  that  they  would  degrade  significantly  over  the  period  of 
a  photographic  run.  In  Figure  IV-1,  runs  were  made  from  low  to  high 
flow  rates  and  back,  as  a  check,  over  which  period  little  degradation 
occurred.  Thus  the  set  of  data  is  not  influenced  by  any  changes  in 
shape  owing  to  degradation. 

In  the  photographic  work  the  data  in  the  figure  were  used  as 
a  reference  for  obtaining  a  consistent  amount  of  degradation.  Data 
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FIGURE  IV  - 1 
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points  obtained  during  the  photographic  runs  are  indicated  by  the  large 
triangles.  Table  I V - 1  summarizes  the  photographic  runs.  Except  for 
run  10  and  12,  it  can  be  observed  that  no  apparent  large  degrada¬ 
tion  occurred  during  the  photographic  runs  from  the  reported  range 
of  friction  velocities  calculated  from  pressure  drop  measurements. 
However,  the  friction  velocities  of  run  12  show  a  significant 
degradation  occurred  during  the  period  of  the  run,  and  consequently  run 
12  was  split  into  two  separate  entities:  run  12a  and  run  12b.  The 
degradation  occurred  simply  because  of  the  inordinate  length  of  time 
necessary  to  complete  the  photographs  for  this  particular  run  at  high 
shear  rate.  For  rur.  10  an  average  friction  velocity  was  used. 

As  noted  from  the  triangles,  drag  reduction  of  up  to  60%  could 
be  obtained  for  the  photographic  runs  in  the  1"  pipe  and  of  44%  in  the 
2.75"  pipe  for  the  dilute  solutions.  The  maximum  of  flow  rate  of  the 
system  prevented  the  photographic  runs  at  higher  N^e  in  the  2.75"  pipe 
while  degradation  caused  by  high  shear  with  vibration  of  the  system 
limited  readings  in  the  1"  pipe. 

The  prediction  of  drag  reduction  in  the  2.75"  pipe  from  the 
1"  pipe  data  points,  using  Equation  11-16  and  the  Seyer  and  Metzner 
B  function,  agreed  well  with  the  experimental  friction  data  as  shown 
on  Figure  IV-1.  These  predicted  lines  were  obtained  from  an  estimated 

i/ 

relaxation  time  function  e  t  obtained  from  the  1"  pipe  friction 

w 

data  (see  Appendix  A).  At  the  highest  Reynolds  number  (highest  shear 
stress)  there  is  a  small  disagreement  between  the  data  points  and  the 
curves,  however,  essentially  exact  agreement  would  have  been  found 
using  point  values  of  e  obtained  from  the  1"  pipe  friction  data. 
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measurements . 

These  calculations  quantitatively  support  the  intercept 
function  determined  by  Seyer  and  Metzner  using  solutions  that  are  sig¬ 
nificantly  more  elastic  and  for  a  moderately  larger  tube  than  they  used. 
In  the  calculations,  however,  the  parameter  e+  (see  Appendix  A) 
is  limited  to  values  less  than  12  so  that  the  calculations  do  not 
provide  a  test  of  the  maximum  drag  reduction  asymptote  as  noted  in 
Chapter  II  (B  =  32).  For  larger  tubes  at  equivalent  Reynolds  number, 
however,  e+  will  generally  be  less  than  12  so  that  predictions  outside 
this  range  are  of  little  consequence. 

In  Figure  I V- 1  ,  the  friction  measurements  of  the  concentrated 
solutions  show  the  characteristic  slow  departure  from  the  laminar 
curve  relatively  to  Newtonian  fluids.  These  solutions  were  prepared 
and  degraded  similarly  to  the  dilute  0.01%  solutions.  The  data  serve 
to  extend  the  previous  measurements  (41)  for  concentrated  drag 
reducers  which  are  known  to  exhibit  transitional  flow  at  surprisingly 
high  Reynolds  numbers.  As  v/ith  the  0.01%,  these  data  show  a  bit  more 
drag  reduction  than  the  same  concentration  studied  by  Seyer  and  Metzner. 
It  is  important  to  notice  that  the  experimental  curves  for  both  the 
1"  and  2.75"  pipes  are  reasonably  close  to  the  predicted  maximum 
drag  reduction  asymptote  of  Virk  et  al  (47)  if  plotted  according  to 

NRe  * 

In  order  to  check  the  asymptotic  value  of  the  B  function, 

the  1"  pipe  data  curve  was  extended,  as  indicated  by  the  dashed  line, 

6 

to  a  Reynolds  number  of  10  ,  where,  as  shown  later,  fully  turbulent 
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flow  is  estimated  to  exist  (estimation  based  on  the  calculated 
intermittency  factor  presented  in  Figure  IV-13).  The  calculated 

r 

value  of  B  at  a  of  10  was  calculated  to  be  43.  This  value  is 
significantly  larger  than  the  asymptotic  value  of  32.2  obtained  from 
Seyer  and  Metzner's  work.  It  should  be  noted  that  the  data  do  not 
prove  the  existence  of  an  asymptote,  but  merely  show  that  previous 
asymptotes  are  not  unique. 
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I V- 2  INSTANTANEOUS  VELOCITIES  AND  ANALYSIS  OF  HISTOGRAMS 

The  data  for  all  runs  of  the  calculated  values  of  velocity 
and  turbulence  intensity  at  a  given  point  in  the  tube  can  be  found  in 
Appendix  H.  These  runs  correspond  to  those  summarized  in  Table  I V- 1 . 

The  instantaneous  velocities  are  available  in  tabulated  form  in 
Appendix  K  enclosed  in  book  II. 

In  order  to  show  that  the  number  of  instantaneous  velocities 
used  to  estimate  the  mean  velocities  was  adequate,  the  cumulative  mean 
of  axial  and  radial  velocities  were  considered  for  each  run.  Figure 
I V - 2  shows,  as  an  example  for  different  radial  positions  in  the  1"  pipe, 
the  levelling  of  the  axial  velocities  for  run  6.  Large  fluctuations 
of  the  cumulative  mean  are  observed  when  less  than  10  velocities  are 
used,  but  the  curves  level  rapidly  for  30  or  more  observations.  It 
is  apparent  from  the  figure  that  with  30  observations  or  more,  random 
fluctuations  in  velocity  have  negligible  influence  on  the  calculated 
mean.  Similar  results  are  predicted  from  simple  statistical  analysis. 

Similarly,  the  typical  cumulative  axial  and  radial  intensi¬ 
ties  (relative  to  the  maximum  velocity)  are  shown  in  Figure  I V - 3 . 

For  these  calculations,  particularly  for  the  axial  intensity,  some  of 
the  curves  tend  to  a  stationary  value  rather  slowly.  This  is  shown 
best  by  the  top  curve  in  Figure  I V- 3 ,  which  does  not  become  stationary 
until  about  70  observations  are  used.  Consideration  of  other  runs 
with  as  many  as  200  observations  show  that  any  change  in  intensity  after 
70  observations  is  insignificant.  Again,  simple  statistics  show  that 
many  more  observations  are  necessary  to  estimate  a  variance  than  a 
mean.  Uncertainty  limits  on  the  statistical  quantities  are  considered 
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FIGURE  IV- 2 
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FIGURE  IV  3 


CUMULATIVE  INTENSITIES 
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in  Appendix  C.  As  illustrated  by  the  figure,  the  axial  intensity  was 
found  to  be  more  sensitive  than  the  radial  intensity. 
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Histograms  of  instantaneous  velocities  for  both  water  and 
0.01%  solution  (run  9  and  run  10),  at  approximately  the  same  Reynolds 
number,  are  given  in  Figures  IV-4  and  IV-5  for  axial  and  radial 
velocities  respectively.  The  gross  features  of  these  histograms 
are  in  qualitative  agreement  with  earlier  independent  studies  and 
offer  support  for  the  bursting  phenomena. 

Except  for  the  histograms  of  radial  velocity  for  polymer, 
the  presence  of  binodal  instantaneous  velocity  distributions,  similar 
to  those  presented  by  Popovich  (33),  can  be  observed  on  Figures  IV-4 
and  IV-5  where  the  histograms  consistently  show  a  dip  near  the  center 
of  the  plot.  These  are  typical  of  all  the  histograms  obtained  in  or 
near  the  boundary  layer.  As  noted  in  Chapter  2,  this  suggests  that 
the  velocities,  inside  the  boundary  layer,  are  composed  of  two  different 
distributions  depending  on  the  origin  of  the  fluid.  The  lower  peak 
can  be  identified  with  the  mean  axial  velocity  of  lumps  of  ejected 
fluid,  while  the  higher  mean  axial  velocity  or  peak  results  from  lumps 
of  fluid  originating  from  the  flow  outside  the  boundary  layer.  Double 
peaks  in  the  core  were  not  observed. 

In  Figure  IV-5,  the  histograms  of  the  radial  velocities  for 
the  0.01%  are  markedly  different  than  the  Newtonian  ones  and  show  no 
indication  of  two  peaks.  This  is  viewed  as  a  consequence  of  the 
resistance  to  stretching  in  the  radial  direction  destroying  the  large 
fluctuations  or  peak  in  the  histograms.  From  the  radial  distri¬ 
butions  it  is  clear  that  the  magnitude  of  the  largest  fluctua- 
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FIGURE  IV- 4 
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FIGURE  IV-  5 
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FIGURE  IV -  6 
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tions  is  reduced  relatively  to  water.  It  is  difficult  to 
make  a  statement  regarding  the  frequency  of  large  fluctuations  and 
whether  they  have  been  reduced  or  not,  as  noted  by  Meek  and  Baer  for 
example  (44),  owing  to  the  small  number  of  large  fluctuations  observed. 

In  order  to  show  in  more  detail  the  nature  of  the  axial 
instantaneous  velocity  distributions,  inside  the  boundary  layer,  the 
distributions  of  Figure  IV-4  have  been  split  according  to  the  sign  of 
the  radial  component  and  are  presented  in  Figure  I V -6  for  run  10. 

The  mean  values  of  axial  instantaneous  velocity  with  positive  radial 
component  were  found  larger  than  the  ones  associated  with  negative 
radial  components.  This  shows  that  the  ejections  must  originate  from 
a  region  of  low  momentum,  and  is  in  agreement  with  Kline  et  al  (20), 
Corino  and  Brodkey  (  '  ),  Klebanoff  (19).  Qualitatively  the  mean  for 
each  side  of  the  distribution  corresponds  to  each  of  the  peaks  in 
the  binodal.  Splitting  of  the  histograms  for  flow  in  core  did  not 
show  that  a  significant  difference  existed  between  means  of  axial 
velocity  associated  with  plus  and  minus  radial  velocities.  The  split¬ 
ting  of  distributions  is  tabulated  in  Appendix  F. 

Comparisons  were  also  made  at  identical  friction  velocity, 
as  suggested  by  previous  work  (15,36,44,16),  and  are  presented  later 
as  intensity  profiles.  In  any  event,  the  radial  fluctuations  of  the 
velocity  were  found  to  be  identically  lower  than  the  Newtonian  fluids 
in  the  boundary  layer  when  normalized  to  the  wall  parameters,  but  the 
same  outside  the  boundary  layer. 

Values  of  the  correlation  factor  R  ,  defined  in  Equation 
1 1 -2 ,  are  presented  as  Figure  I V- 7  for  runs  6  and  10  of  the  dilute 
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polymer  and  run  9  for  water.  Data  for  the  other  runs,  which  have  been 
omitted  to  clarify  the  figure,  show  similar  behavior.  Although  there 
is  considerable  scatter  in  the  data,  there  are  distinct  trends  with 
radial  position.  As  discussed  earlier,  the  bursting  arguments  coupled 
with  simple  continuity  considerations,  imply  Rq  should  be  greater  than 
unity  in  the  vicinity  of  the  edge  of  the  sub-layer.  The  data  points 
for  the  water  connected  with  the  dashed  line,  although  not  extending 
deeply  enough  into  the  sublayer,  suggest  this  trend.  To  map  the 
trend  completely,  one  would  like  measurements  for  Y+  <  10  or  to  posi¬ 
tions,  as  indicated  previously  (/  ),  where  the  ejections  originate. 

For  example  for  water,  the  data  point  closest  to  the  wall,  at  a 
Y+  -  14,  has  Rq  of  approximately  3.  Equivalently,  at  this  position 
positive  radial  fluctuations  are  observed  three  times  as  often  as  the 
larger  negative  fluctuations,  while  the  mean  of  all  the  observed 
fluctuations  is  zero  within  the  statistical  uncertainty  of  the 
cal cul ations . 

For  the  polymer  solutions,  runs  6  and  10,  there  is  a  distinct 
peak  in  values  of  Rq  at  £  of  about  0.05.  The  values  become  constant 
at  about  £  =  0.1.  Consideration  of  the  dimensionless  velocity  pro¬ 
files  shown  later  in  Figure  I V- 1 6 ,  indicates,  as  expected,  this  peak 
is  within  the  sublayer  and  the  end  of  the  peak  (c-0.1)  coincides  with 
the  outer  edge  of  the  sublayer  which  has  been  shifted  to  about  Y+  -  100 
for  the  polymer.  Qualitatively  the  bursting  arguments  suggest  that 
R  ,  in  addition  to  being  greater  than  unity,  should  show  a  peak  owing 
to  the  acceleration  from  small  to  large  radial  velocity  as  the  element 
of  fluid  moves  through  the  sublayer.  The  values  should  decrease  to 
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unity  at  the  edge  of  the  sublayer  as  the  fluid  mixes  with  the  core 
fl  uid. 

For  the  remainder  of  the  cross  section,  for  all  the  runs, 
Rq  scatters  around  unity,  except  over  a  narrow  range  of  radial  posi¬ 
tions  near  to  5  -  0,65  where  it  is  consistently  less  than  unity, 
Seyer's  earlier  data  (41)  show  exactly  the  same  trends,  however,  no 
explanation  can  be  offered  for  this  behavior. 
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IV— 3  VELOCITY  PROFILES 

Figures  IV-8  to  IV-10  represent  actual  time  average  velocity 
profiles  for  all  the  runs  in  the  1"  and  2.75"  pipes.  The  mean  or  bulk 
velocities  obtained  by  integrating  these  profiles  (using  actual  data 
points  and  not  smoothing  the  profiles)  agreed  within  2%  of  the  recorded 
bulk  velocities,  except  for  run  12  which  was  4%  in  error  (see  Table 
I V- 1 ) .  Owing  to  the  small  amount  of  scatter,  these  data  show  the  use¬ 
fulness  of  the  streak  photograph  method  over  a  wide  range  of  conditions. 
It  should  be  noted  that  no  discontinuity  appears  in  the  curves  owing 
to  different  magnifications  used  during  a  run. 

More  information  about  the  nature  of  the  profiles  is  found 
in  Figure  IV-11  where  dimensionless  velocity  profiles  are  compared  to 
predictions  using  Bogue's  equation  for  turbulent  flow  and  the  Newtonian 
laminar  equation.  The  water  data  are  in  excellent  agreement  with  the 
predicted  curve  while  the  0.01%  profile  which  is  at  the  highest 
Reynolds  number  for  the  dilute  solutions  is  markedly  flatter  in  the 
core.  For  the  other  runs  at  lower  Reynolds  number  where  less  drag 
reduction  is  occurri ng,  the  flattening  is  not  so  large  and  the  curves 
approach  Newtonian  behavior. 
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FIGURE  IV-  8 
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FIGURE  IV- 9 
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FIGURE  IV-10 
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FIGURE  IV-11 
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IV-4  TRANSITIONAL  FLOW 

Figure  I V- 1 2  is  a  plot  of  the  velocity  data  for  the  0.10% 
solutions  which  serve  to  extend  Seyer's  earlier  data  (41  )  for  transi¬ 
tional  flow.  The  fact  that  these  profiles  reflect  transitional  flow 
has  been  documented  in  detail  before.  Briefly,  the  transitional  flow 
consists  of  alternating  patches  of  turbulent  and  laminar  fluid.  The 
relative  number  of  each  is  described  by  the  i ntermi ttency  factor,  y, 
defined  as  the  fraction  of  the  time  the  flow  is  turbulent.  Thus  y 
equal  unity  is  fully  turbulent  while  y  of  zero  reflects  laminar  flow. 
Furthermore,  it  has  been  shown  (35)  that  the  i ntermi ttency ,  at  the 
center  of  the  pipe,  is  given  simply  by  the  position  of  the  velocity 
profile  relative  to  the  fully  turbulent  and  fully  laminar  profiles. 

Figure  I V- 1 3  represents  the  i ntermi ttency  factors  calculated 
from  the  velocity  profiles  at  the  center  line  of  both  1"  and  2.75" 
pipes  and  are  compared  with  Rotta's  data  points  for  water  (35).  The 
1"  pipe  intermi ttency  factor  curve  indicates  how  slowly  y  increases 
with  the  Reynolds  number  for  the  concentrated  drag  reducer  as  compared 
to  the  Newtonian  curve  (see  Rotta  (35  )  and  Patel  and  Head  (31  )). 

This  curve  indicates  also  that  the  fully  turbulent  flow  will  be  reached 

r 

at  a  Nj^e  -  10°.  Although  there  is  some  uncertainty  in  locating  the 
point  of  zero  i ntermi ttency ,  it  is  clear  that  the  Reynolds  number  where 
turbulence  first  appears  is  not  significantly  different  from  2100. 

From  the  frictional  data  curve  on  Figure  IV-1 ,  it  was  not  so  evident 
that  the  transitional  flow  existed  at  that 

The  velocity  measurements  indicated  that  the  intermi ttency 
factor  depends  strongly  on  the  tube  diameter  where  a  39%  increase 
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FIGURE  IV-13 
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relative  to  the  1"  pipe  curve  was  recorded  for  the  2.75"  pipe  at  a 

4 

NRe  of  2  x  10  .  Similarly  the  i ntermi ttency  factors  determined  from 
Seyer's  profiles  in  a  1"  pipe  are  significantly  larger  than  the  ones 

5 

in  this  work  and  suggest  a  Reynolds  number  of  about  10  is  necessary 
for  fully  turbulent  flow.  Thus  the  i ntermi ttency ,  at  high  Reynolds 
number  is  very  sensitive  to  the  level  of  degradation  of  the  polymer. 

Figures  I V- 1 4  to  IV-16  give  the  corresponding  intensities 
for  the  transitional  flow.  These  simply  confirm  Seyer's  earlier 
data  (41)  which  shows  a  very  small  radial  intensity  but  an  apparently 
large  axial  intensity.  The  radial  fluctuations  decrease  rapidly  to 
zero,  because  they  can  only  occur  y  of  the  time;  however,  large  axial 
fluctuations  still  occur  because  of  the  profile  alternating  between 
fully  laminar  and  fully  turbulent  profiles. 

These  measurements  make  it  apparent  that  the  limits  of 
transitional  flow  should  be  investigated  in  order  to  determine  the 
range  of  applicability  of  the  similarity  laws  for  drag  reducers.  The 
limits  must  be  related  to  the  thickening  of  the  boundary  layer  as 
defined  by  Equation  11-12  for  example,  and  clearly  must  be  less  than 
the  tube  radius.  Currently  the  true  upper  limit  of  stable  boundary 
layer  thickness  is  unknown;  however,  for  Newtonian  fluids  at  Reynolds 
number  of  2100  the  thickness  is  approximately  £  =  0.1. 
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IV— 5  SEMI  LOGARITHMIC  VELOCITY  PROFILES 

Velocity  profiles  for  the  dilute  drag  reducer  and  water, 

plotted  on  U+  -  Y+  coordinates  are  given  by  Figure  IV-17.  For  the 

polymer  solution  the  solid  lines  represent  eyeball  estimates  of  the 

best  curve  through  the  data.  The  curve  through  the  water  data  has 

been  determined  using  Equation  11-11.  At  Reynolds  number  of  80354 

in  the  1"  tube  there  is  a  considerable  shift  in  the  data  points  for 

30  <  Y+  <  ioo.  The  upper  points  which  form  a  continuous  curve  to 

Y+  about  90  were  obtained  at  high  magnification  while  the  remaining 

points  were  obtained  with  a  lower  magnification.  As  noted  before 

and  in  Table  IV-1  with  runs  designated  12a  and  12b  this  run  suffered 

from  a  significant  amount  of  degradation.  Thus  the  shift  simply 

* 

reflects  the  increase  in  U  that  occurred  between  the  high  and  low 

magnification  portions  of  the  run.  In  the  other  runs  where  a  change 
★ 

in  U  was  noted  the  average  was  used  in  Figure  IV-17  and  one  cannot 
distinguish  between  data  obtained  at  high  and  low  magnifications. 

In  all  cases  the  data  correspond  roughly  to  Equation  1-1 
but  not  in  detail.  In  particular  the  Newtonian  data  have  a  constant 
intercept  B  and  slope  A  up  to  Y+  about  400  but  deviate  considerably 
at  higher  Y+.  This  behavior  is  quantitatively  equivalent  to  that 
observed  in  the  number  of  other  studies  reporting  turbulent  velocity 
profiles  of  Newtonian  fluids  (18). 

The  velocity  profiles  for  the  0,01%  polymer  solutions  show 
a  marked  similarity  to  the  Newtonian  data  at  high  Y+  suggesting  the 
same  correction.. Cg^f ) ,  may  be  applicable.  From  the  figure  it  is 
evident  that  the  velocity  profiles  could  be  interpreted  in  a  number  of 
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FIGURE  IV-17 
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ways : 

(i)  using  Equation  1-1  the  best  fit  line  through  the  data  would 
clearly  have  a  slope  somewhat  higher  than  the  value  2.46  for 
Newtonian  fluids. 

(ii)  using  Equation  1-1  with  A  =  2.46  would  fit  the  data  poorly, 
at  least  at  the  higher  Reynolds  numbers,  and  would  result  in 
a  large  uncertainty  in  vertical  location  of  the  line. 

(iii)  using  Equation  1 1 - 1 1  and  assuming  Bogue's  correction  function, 
at  least  as  a  first  approximation,  is  applicable  to  the  drag 
reducers . 

It  is  emphasized  that  each  of  the  above  yield  different 
values  of  the  intercept  function  B.  In  the  extreme  case  the  data  in 
Figure  IV-17  shows  that  the  discrepancy  between  (i)  and  (iii)  above 
could  be  as  large  as  20%.  In  the  case  of  predictions  for  small 
amounts  of  drag  reduction  in  large  tubes  an  error  of  20%  in  B  can  be 
shown  to  cause  significant  discrepancies  in  predicted  pressure  drops. 

For  drag  reducers  apparently  only  one  study  has  applied 
Bogue's  correction  function  to  velocity  profiles.  As  mentioned 
earlier  Ni codemo  et  al  (29)  have  measured  velocity  profiles  in  typical 
drag  reducers  using  a  Pitot  probe  calibrated  in  a  circular  tow  tank 
filled  with  the  polymer  solution  under  consideration.  Based  on  the 
observation  that  the  data  on  a  semilog  plot  are  linear  over  only  a 
small  portion  of  the  tube  near  the  center,  they  conclude  the  profiles 
cannot  be  interpreted  according  to  Equation  11-11  without  major 
modification.  In  contrast  the  velocity  profile  obtained  by  Seyer  (41) 
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from  streak  photographs  is  in  good  agreement  with  Equation  11-11.  A 
detailed  examination  of  the  profile  however,  indicates  the  data  are 
really  inconclusive  and  could  be  just  as  well  interpreted  with  or 
without  the  function  C^U, f).  This  being  mostly  because  Cg(s,f)  is 
negligible  for  Seyer's  data. 

In  Figure  IV-17  Bogue's  correction  function  has  been  applied 
to  the  smoothing  curves  drawn  through  the  polymer  data.  The  resulting 
corrected  velocity  profiles  are  given  by  the  lines  located  below  the 
data  points  of  the  various  runs.  For  all  of  the  runs  the  correction 
applied  to  the  smoothed  curve  yields  a  line  that  is  parallel  (within 
the  uncertainty  of  locating  the  smoothed  curve)  to  the  corrected 
Newtonian  curve  showing  that  (a)  Bogue's  correction  function  is 
adequate  for  the  drag  reducers,  and  (b)  the  constant  A  is  the  same  as 
for  Newtonian  fluids. 

Examination  of  the  resulting  velocity  profiles  shows  that 
the  difficulty  in  interpreting  the  uncorrected  profiles  results  from 
the  greatly  thickened  transition  portion  of  the  profile  merging  with 
the  high  Y  data  approximately  where  the  correction  function  becomes 
significant.  For  Newtonian  fluids  there  is  a  large  range  of  Y+(30  < 

Y+  <  400)  separating  these  two  regions. 

Figure  IV-17  shows  quantitatively  the  marked  thickening  of 
the  sublayer  region.  For  the  highest  Reynolds  number  the  influence 
of  the  wall  extends  to  Y  =  200  whereas  for  Newtonian  fluids  Y  =  30 
is  about  the  upper  limit  of  the  wall  region.  It  is  interesting  to 
note  that  the  buffer  region  is  somewhat  larger  in  the  2.75"  than 
that  in  the  1"  tube. 


. 


If  the  correction  factor  is  ignored  significant  errors  may 
result  in  the  determination  of  the  intercept  B(e+)  and  uncertainty  in 
the  intercept  function  may  be  serious  in  prediction  of  drag  reduction 
particularly  in  large  pipes. 

Since  no  elasticity  measurements  were  made  in  this  work  it 

is  not  possible  to  determine  B  as  a  function  of  the  parameter  e+ 

as  is  suggested  by  Seyer  (43).  However,  for  a  given  polymer  and 

+  * 

concentration  0  depends  on  U  and  therefore  it  is  equivalent  to  plot 

•k 

B  as  a  function  of  U  . 

The  function  B  as  determined  from  friction  measurements  and 
corrected  velocity  profiles  is  given  by  Figure  IV-18.  The  agreement 
from  the  velocity  profiles  of  the  two  tubes  as  well  as  that  from 
friction  measurements  supports  in  some  detail,  and  over  a  significant 
range  of  variables,  the  velocity  profile  and  attendent  analysis 
given  by  Equations  II-l  to  1 1 - 1 1 .  The  largest  disagreement  occurs  for 
the  point  designated  as  run  10.  As  indicated  by  the  somewhat  lower 
friction  factor  for  this  point  in  Figure  IV-1 ,  the  polymer  had  not 
been  degraded  sufficiently  and  so  the  higher  value  of  B  is  expected. 
Alternatively,  if  the  function  were  plotted  against  the  relaxation 
time  parameter  the  higher  elasticity  of  the  material  at  conditions 
of  run  10  would  cause  this  point  to  shift  to  the  right  relative  to 
points  6  and  16,  and  would  fall  much  nearer  to  the  curve. 
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FIGURE  IV-18 
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Rudd  (37)  indicated  that  if  Schubert  and  Corcos  theory  of 
the  wall  turbulence  is  considered  and  modified  for  the  case  of  visco¬ 
elastic  fluid,  then  the  similarity  laws  can  be  rescaled  to  take 
account  of  the  viscoelastic  terms  so  that  the  parameters  for  drag 
reducers  can  be  predicted  from  the  Newtonian  values.  From  his  develop¬ 
ment  and  experimental  determination  of  the  scale  factor,  he  obtained 
the  following  expression  for  the  B  function 


B  =  12  D  -  2.44  1 n D  -  6.5  IV-1 

where 

D2  =  1  +  K  e+  I V - 2 


is  the  scale  factor  with  K  -  0.5  as  found  from  Rudd's  data.  This  scale 
factor  should,  if  calculated  for  each  run,  rescale  all  turbulence 
parameters  (velocity  and  distance)  to  bring  them  on  a  single  curve, 
the  Newtonian  case.  From  his  work  the  predicted  values  of  B  with 
e+  were  compared  to  the  experimental  results  of  Seyer  and  Metzner. 

The  B  functions  are  in  agreement  for  e+  lower  than  17  but  Rudd's  esti¬ 
mated  asymptote  is  higher  with  value  of  40  for  e+  of  37.  This  is  in 
reasonable  agreement  with  the  extrapolation  of  the  friction  data  in 
this  work. 

It  is  interesting  to  use  his  correlation  to  estimate  the 
boundary  layer  thickness  and  the  value  of  B  for  run  12b  of  this  work. 
The  estimated  relaxation  time  obtained  from  the  friction  data  of  run 

12b  and  using  Equation  11-16  was  used  in  Equation  IV -2  to  evaluate 

-4 

the  scale  factor  D.  It  was  found  that  0  =  8.8  x  10  sec  and  D  =  2.4. 
The  estimated  dimensionless  thickness  of  the  boundary  was  Y+  =  84 
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(Y,  x  D)  and  the  value  of  B  was  found  to  be  20.61.  From  Figure 
water 

IV-17  the  intercept  is  19.1  for  this  run  and  the  boundary  layer 
extends  to  something  in  the  order  of  100.  This  is  surprisingly 
good  agreement. 
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I V- 6  TURBULENCE  INTENSITIES 

The  RMS  values  obtained  for  each  set  of  the  axial  and  radial 
instantaneous  velocities  are  used  to  estimate  the  relative  turbulence 
intensities.  As  indicated  in  Chapter  I,  great  confusion  is  found  from 
a  survey  of  the  measured  relative  intensities  in  the  drag  reducers 
owing  to  the  uncertainty  in  accuracy  of  the  results  and  in  how  to 
scale  or  compare  them  for  different  systems. 

The  sample  histograms  for  polymer  in  Figure  I V -5  when  compared 
to  water  have  already  indicated  a  significant  lowering  of  radial 
fluctuations  or  intensity  occurs  for  the  dilute  drag  reducer  at  a  given 
Y+. 

Figures  IV-1  9  to  IV-2  3  show  all  of  the  Newtonian  intensity 
data  obtained  in  this  work.  The  curves  on  the  figures  represent 
Sandborn's  data  for  air  at  the  indicated  Reynolds  numbers.  Runs  4  and 
5  on  Figures  IV-1 9  and  IV- 20  do  not  extend  over  an  appreciable  portion 
of  the  tube  but  are  considered  to  be  in  excellent  agreement  with  the 
curves  at  the  smallest  radial  positions.  Reference  to  Appendix  H  shows 
these  data  are  based  on  a  small  number  of  instantaneous  readings 
relative  to  the  other  runs  in  this  work. 

In  Figure  IV-21  for  run  9  the  circled  points  form  a  group  of 
data  which  are  consistently  lower  than  the  curves.  These  data  which 
have  been  circled  are  from  a  separate  set  of  photographs  (the  high 
magnification  portion  of  the  run)  than  the  remaining  points  and 
apparently  reflect  some  mistake  in  determining  calibration  constants. 

If  they  are  ignored  the  data  show  excellent  agreement  with  Sandborn's 
curves  over  the  entire  section. 
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The  data  in  Figures  IV-22  and  IV-23  for  the  remaining  runs 
which  are  at  the  highest  Reynolds  numbers  are  again  in  excellent  agree¬ 
ment  with  the  curves  over  most  of  the  section.  Near  the  centerline 
the  data  are  somewhat  higher  than  the  curves,  however,  it  is  not  clear 
what  the  origin  of  this  discrepancy  is.  As  shown  in  Appendix  E, 
calculations  of  the  Reynolds  stress  by  forming  the  product  U'U'  and  the 

A  I 

kinetic  energy  U'U1  for  run  11  are  in  excellent  agreement  with  the 
measured  Reynolds  stress  and  Laufer's  (21)  kinetic  energy.  Thus  it  is 
not  clear  that  the  intensities  are  seriously  in  error. 
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The  .01%  intensities  relative  to  maximum  velocity  for  runs 
6,  10  and  12  in  the  1"  tube  and  run  16  in  the  2.75"  tube  are  compared 
to  the  air  curves  in  Figures  IV-24  to  IV-27.  These  figures  provide  a 
comparison  of  the  intensity  calculated  relative  to  maximum  velocity  at 
approximately  equal  Reynolds  number.  Since  the  water  and  polymer 
viscosities  are  about  the  same  the  comparison  could  also  be  viewed  as 
one  at  the  same  flow  rate  or  bulk  average  velocity. 

In  each  run  the  radial  intensities  are  markedly  lower  at 
all  radial  positions  than  the  intensities  for  Newtonian  fluids.  The 
amount  of  lowering  of  the  intensities  is  ordered  according  to  the 
amount  of  drag  reduction;  thus  the  radial  intensities  for  the  2.75" 
pipe  are  somewhat  higher  than  those  in  the  1"  pipe  at  an  equivalent 
Reynolds  number. 

Axial  intensities  for  the  low  Reynolds  number  runs, 
run  6  and  run  10  are  not  altered  significantly  from  the  Newtonian 
values.  At.  high.  Reynolds  number  however,  they  are  significantly 
lower  in  the  core.  Notice  that  the  Newtonian  curves  of  axial  intensity 
in  Figures  IV-26  and  IV-27  are  for  Reynolds  number  of  170000  and  that 
a  curve  for  Reynolds  number  coinciding  with  the  data  points  would  be 
above  this  curve. 

In  all  cases  the  rapid  increase  in  intensity,  usually 
associated  with  the  edge  of  the  boundary  layer,  is  shifted  toward 
larger  radial  positions.  This  suggests,  as  do  the  velocity  profiles, 
a  thickening  of  the  boundary  layer  occurs. 
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The  intensity  data  made  dimensionless  with  friction  velocity 
*  + 

U  are  plotted  as  a  function  of  Y  in  Figure  IV-28.  The  effect  of 

k 

dividing  by  U  is  to  shift  the  intensities  of  the  drag  reducers  upward 
relative  to  the  water  curves.  In  the  lower  portion  of  the  figure  the 
radial  intensities  still  tend  to  be  lower  than  the  Newtonian  intensi¬ 
ties,  particularly  in  the  wall  region.  The  general  trend  is  that  the 
polymer  intensities  are  increasing  more  slowly  with  radial  position, 
in  accord  with  the  thickening  of  the  sublayers.  Ultimately,  they 
attain  about  the  same  value  as  for  Newtonian  fluids.  In  terms  of  the 
bursting  arguments  the  radial  intensities  indicate  that  the  large  wall 
eddies  are  simply  rescaled  in  size.  It  is  apparent,  however,  since 
the  frequency  of  occurrence  can  be  reduced,  as  shown  in  studies  by 
Fortuna  and  Hanratty  (15)  as  well  as  Meek  and  Baer  (/a),  that  the 
process  of  lowering  the  radial  intensities  is  more  complicated  than 
scaling  just  the  size. 

In  contrast  to  the  radial  intensities  the  axial  intensities 

k 

in  Figure  IV-28  relative  to  U  are  all  significantly  larger  than  the 
Newtonian  values  in  the  wall  region  and  approach  a  single  curve  in  the 
core  region.  The  shape  of  these  curves  is  identical  to  Rudd's  (36) 
for  a  similar  amount  of  drag  reduction.  However,  these  data  are  not 
ordered  with  respect  to  amount  of  drag  reduction  on  this  type  of  plot. 
Reference  to  Table  IV-2  shows  that  runs  5,  6,  10  in  the  1"  tube  and 

k 

16  in  the  2.75"  tube  are  all  at  approximately  the  same  U  (the  value 
of  U  for  run  12  is  two  times  larger). 

In  Figures  IV-17  and  IV-18,  it  is  seen  that  these  runs 
have  increasing  intercepts  owing  to  different  levels  of  elasticity 
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of  the  fluid.  The  increased  intensities  are  in  the  order  of  increas¬ 
ing  elasticity  or  equivalently  increase  in  the  velocity  profile 

•k 

intercept,  B,  at  fixed  U  .  It  is  obvious  that  the  level  of  intensity, 
* 

at  a  fixed  U  ,  is  very  sensitive  to  the  elasticity  of  the  system. 

Consideration  of  the  velocity  intercept  for  run  12,  which 
has  the  highest  intercept  of  all  the  runs,  would  indicate  the  inten¬ 
sities  for  this  run  should  be  extremely  high.  This  is  not  the  case, 
however,  as  these  intensities  are  lower  than  for  runs  10  and  16  and 
about  the  same  as  those  of  run  6.  In  the  core  these  data  are  con¬ 
sistent  with  the  other  runs. 

Dimensional  considerations  suggest  (see  Chapter  II)  that 
the  increase  in  intensities  in  the  wall  region  should  be  ordered 

'kO 

+  0  u  ^ 

with  respect  to  the  parameter  e  =  — —  .  Rudd  (3/)  as  indicated 

earlier  has  formulated  an  explicit  expression  for  the  dependence 

of  the  flow  parameters  on  e+.  Figure  IV-29  is  a  plot  of  the  rescaled 

intensities  and  distances  as  suggested  by  Rudd.  The  parameter  D, 

defined  in  Equation  IV-1  ,  has  been  determined  directly  by  substituting 

values  of  B  from  the  velocity  profiles.  Notice  if  equation  IV-2  were 

used  to  find  relaxation  times,  the  times  so  determined  would  be 

slightly  different  than  those  shown  in  Appendix  A.  The  axial  intensi- 

★ 

ties  for  the  runs  with  equal  U  are  collapsed  to  single  curve 
which  coincides  with  the  curve  for  water.  Rudd  has  based  his  con¬ 
clusions  on  a  similar  plot  for  only  one  intensity  profile.  The 

* 

intensities  for  run  12,  which  is  at  different  U  ,  are  not  scaled 
properly.  In  fact  the  spread  of  intensities  between  run  12  and  the 
master  curve  is  about  the  same  as  the  maximum  spread  in  the  previous 
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The  radial  intensities  do  not  superimpose  on  the  water 
curve,  but  are  much  lower  at  all  radial  positions.  The  polymer  inten¬ 
sities  follow  a  single  curve  but  there  is  no  obvious  improvement  over 
the  plot  of  the  same  data  in  Figure  IV-28. 
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CONCLUSION  AND  RECOMMENDATIONS 


An  extensive  series  of  measurements  of  instantaneous  velocity 
streak  photographs  have  been  obtained  in  Newtonian  and  drag  reducing 
fluids.  Based  on  a  comparison  with  available  Newtonian  data  of  velocity 
profiles  and  axial  and  radial  intensities,  the  quantitative  usefulness 
of  the  technique  is  verified.  Since  the  streak  photograph  technique 
does  not  suffer  from  the  serious  limitations  of  probe  devices  in  drag 
reducers  it  is  implicit  that  the  measurements  in  the  drag  reducing 
fluids  are  also  quantitatively  correct. 

Histograms  of  instantaneous  axial  velocity  confirm  the 
presence  of  binodal  distributions  in  the  boundary  layer.  Histograms 
of  radial  velocity  are  binodal  for  water  but  there  is  no  evidence 
of  two  peaks  for  the  polymer  solutions.  By  ordering  the  set  of 
instantaneous  velocities  according  to  the  sign  of  radial  component 
it  is  found  that  on  the  average  the  fluctuating  velocities  toward 
the  center  of  tube  are  large  than  those  toward  the  wall.  This 
observation  agrees  with  the  visual  observations  of  "bursting"  described 
by  Brodkey  amongst  others.  For  a  given  radial  position  the  largest 
fluctuations  in  radial  histograms  for  polymer  solutions  are  less  than 
those  for  Newtonian  fluids  at  similar  conditions. 

For  drag  reducing  fluids  Bogue's  correction  function  must 
be  aoplied  to  the  data  in  order  to  obtain  a  consistent  representation 
of  velocity  profiles.  Failure  to  use  the  correction  function  can 
result  in  large  uncertainties  in  det.ermininq  the  intercept  B(e+)  of 
the  velocity  profiles. 


' 
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The  correct  form  of  the  function  B(e+)  determined  by  Seyer 
and  Metzner  is  correct  for  e+  less  than  approximately  12.  At  large  e+, 
where  a  drag  reduction  asymptote  is  predicted  the  function  is  not  unique. 
The  data  in  this  work  show  that  lower  friction  factors  can  be  observed 
than  Seyer  and  Metzner's. 

Turbulence  i ntermi ttency  factors,  determined  from  velocity 
profiles  for  a  0.1%  solution  in  a  one-inch  pipe  show  the  first 
turbulence  appears  at  a  generalized  Reynolds  number  of  2100.  Extra¬ 
polation  of  the  data  show  that  the  flow  is  not  fully  turbulent  until 
Reynolds  numbers  in  the  vicinity  of  10^  are  attained.  Consideration 
of  earlier  data  and  data  obtained  in  a  2.75"  pipe  show  the  i ntermi ttency 
is  highly  sensitive  to  level  of  elasticity  as  well  as  the  diameter. 

These  measurements  show  the  need  for  further  consideration  of  the  so- 
called  asymptote  as  well  as  the  limits  of  transition. 

If  the  root  mean  square  of  the  axial  and  radial  velocity 
fluctuations  of  the  polymeric  solutions  are  made  dimensionless  with 
the  maximum  velocity,  both  are  found  to  be  lower  than  those  for  water 
at  all  radial  positions;  the  amount  of  lowering  being  larger  as  the 
level  of  elasticity  of  the  solutions  increases.  On  the  other  hand  when 
compared  to  friction  velocity  the  axial  fluctuations  are  larger  while 
the  radial  remain  smaller  than  those  for  Newtonian  fluids.  Scaling  of 
the  intensities  according  to  Rudd  does  not  superimpose  the  curves  unless 
only  data  obtained  at  the  same  friction  velocity  are  considered. 

Further  work  is  needed  to  determine  the  proper  scaling  factors 
for  turbulence  in  these  systems. 


■ 
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RHEOLOGICAL  BEHAVIOR  OF  THE  SOLUTIONS 
i )  0.01%  Solution 

A  sample  of  solution  was  drawn  from  the  experimental  apparatus 
when  pressure. drop  measurements,  recorded  in  the  1"  and  2.75"  pipes, 
corresponded  to  the  predetermined  curves  for  0.01%  solution  as  presen¬ 
ted  in  Figure  IV-l .  Viscometric  measurements  were  done  with  a 
Weissenberg  Rheogoniometer  precalibrated  with  a  standard  oil.  Curves 
for  the  standard  oil  and  polymer  solutions  are  presented  in  Figure  A-l. 
The  straight  line  obtained  for  the  standard  oil  is  representative  of 
Newtonian  behavior  with  a  slope  of  unity  and  a  viscosity  of  8.62  cps. 

The  specified  viscosity  at  the  same  temperature  is  8.63  cps. 

The  curve  of  the  0.01%  solution  indicated  a  slightly  pseudo¬ 
plastic  behavior  in  the  range  of  shear  rates  obtained  with  the  Rheogoni¬ 
ometer.  The  flow  behavior  index  varied  from  0.922  to  0.930  and  the 
apparent  viscosity  from  1.25  cps  to  1.50  cps  at  the  lowest  shear  rate. 
Although  the  Rheogoniometer  data  are  at  much  lower  shear  rates  than 
those  in  the  pipe  experiments,  they  join  smoothly  with  Seyer's  data  (41) 
at  the  shear  stresses  of  interest.  Thus  the  kinematic  viscosity  is 
assumed  to  be  1.2  x  10"^  ft^/sec  for  tw  ranging  from  0.06  to  0.25  lb/ft^, 
which  is  the  value  used  by  Seyer.  Table  A-l  represents  the  viscometric 
data  and  the  results  obtained  from  the  Rheogoniometer  for  the  standard 


oil  and  0.01%  solution. 


' 
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FLOW  CURVES  FOR  POLYMERIC  SOLUTIONS 


TABLE  A- I 


VISCOMETRIC  MEASUREMENTS 
FOR  OIL  AND  POLYMERIC  SOLUTIONS 


Standard  Oil  (9  centipoise)  at  22.5°C 


RUN 

DISPLACEMENT 

ROTATION 

SHEAR  RATE 
sec  1 

SHEAR 

STRESS 

volts 

rpm 

dynes/cm2 

psf 

1 

8.85 

5.68 

35.80 

3.04 

0.00637 

2 

282.00 

17.09 

113.20 

9.68 

0.02030 

3 

448.00 

28.40 

179.00 

15.38 

0.03220 

4 

44.50 

2.84 

17.90 

1.53 

0.00321 

5 

14.00 

0.90 

5.68 

0.48 

0.00101 

Separan  AP 

-30  0.01%  at 

c_> 

o 

LO 

CXI 

CXI 

1 

9.00 

2.84 

17.90 

0.31 

0.00065 

2 

25.00 

9.00 

56.80 

0.857 

0.00179 

3 

72,00 

28.40 

179.00 

2.470 

0.00515 

4 

111.00 

45.20 

285.00 

3.800 

0.00793 

5 

38.60 

14.30 

90.00 

1.300 

0.00271 

6 

71  .00 

28.40 

179.00 

2.440 

0.00509 

RHEOGONIOMETER  PARAMETERS 

a 

-  angle  of  the  cone 

0.952° 

d 

-  diameter  of  the  cone 

10.00  cm 

Kt 

-  torsion  bar  constant 

8.916  dynes-cm/microns 

At 

-  displacement 

volts 

3 

-  transducer  calibration 

1 .009  micron/vol t 

W 

-  rotational  speed 

rpm 

FORMULAE 

SHEAR  STRESS  x  =  (3.82  3  At  Kt )/d3  =  0.0343  At  dynes/cm2 

SHEAR  RATE  8<u"  >/D  =  (360  W)/ (60  a)  =  3.61  W  sec'1 

A 


■ 

• 
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i i )  0.10%  Solution 


The  experimental  points,  calculated  from  the  laminar  friction 
data,  are  compared  to  Seyer's  (41)  viscometric  data  for  the  0.10%  solu¬ 
tions  on  Figure  A-l.  Calculated  values  are  tabulated  in  Table  A- 1 1 . 

The  procedural  similarity  in  the  preparation  of  the  solutions  and  the 
excellent  agreement  obtained  with  Seyer's  data  for  the  same  concentra¬ 
tion  show  the  shear  rate-shear  stress  behavior  of  the  0.10%  solution  is 
the  same.  Therefore  for  cal culational  purposes  the  same  relation 
assumed  by  Seyer  will  be  used  for  6  x  10  <  ^  <  1.2  x  10  . 

8<u>  8<u  >  o 

In  tw  =  -6.540  +  0.598  In  — +  0.00599  (In  — ^-)c 

where  x, ,  is  in  psf 
w  r 


8<v 

D 


.  .  -1 
is  in  sec  . 


TABLE  A- 1 1 


EXPERIMENTAL  VISCOMETRIC 

Flow  (US  gal/min) 

Manometer  Reading  (inches  Hg) 
Length  Between  Pressure  Taps  (ft) 
Diameter  of  Pipe  (ft) 

Temperature  (°F) 

Bulk  Velocity  (ft/sec) 

Pressure  Drop  (psf) 

Shear  Stress  at  Wall  (psf) 

Shear  Rate  ( 8<lT  >/D)(sec  ^) 

X 


DATA  FOR  0.10%  SOLUTION 


5.0 

11.0 

26.4 

27.5 

1.31 

2.09 

3.42 

3.68 

17.5 

17.5 

17  5 

17.5 

0.230 

0.230 

0.230 

0.230 

74  0 

74.0 

74,0 

74.0 

0.267 

0.591 

1 .419 

1  .478 

2.073 

3.308 

5,413 

5.824 

0.0060 

0,0109 

0.0178 

0.0191 

9.32 

20.56 

49,36 

51 .43 
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i i i )  Estimate  of  the  Elasticity  Level  and  Calculations  for 

Prediction  of  2.75"  Friction  Factor  from  1"  Data 


Using  Equations  11-16,  11-12  and  1-4,  four  of  the  experimental 
friction  points  tabulated  in  Appendix  I  for  the  0.01%  solution  were 
used  to  estimate  (by  a  trial  and  error  solution  of  Equation  11-16) 
values  of  B  and  consequently  relaxation  time  in  the  range  of  the  experi¬ 
mental  photographic  runs.  These  calculations  are  tabulated  in  Table 
A- III.  Figure  A-2  shows  the  relaxation  time  as  a  function  of  tw  along 
with  Seyer's  (41)  and  Oliver's  (30)  data  for  their  respective  0.01% 
Separan  solutions.  As  suggested,  from  the  lower  friction  measurements, 
the  level  of  elasticity  of  the  present  solutions  is  higher  than  Seyer's 
solutions.  The  dashed  line  is  taken  to  represent  the  relaxation  time 
of  the  0.01%  solution  in  this  work  and  for  purposes  of  calculation, 
can  be  represented  empirically  by 

In  6  =  0.000388  -  0.4930  In  t 

w 

where  tw  is  expressed  in  psf  and  e  in  sec  \ 

TABLE  A- II I 


LEVEL  OF  ELASTICITY  FOR  0.01%  SOLUTIONS 


Reynolds  number 

132980 

93369 

50000 

30000 

Friction  factor 

0.00140 

0.00177 

0.00270 

0.00382 

B(e+) 

23.96 

20.30 

15.00 

11.50 

+ 

e 

11.70 

9.25 

6.00 

4.00 

Shear  stress  at 

wall  (psf)  0.496 

0.309 

0.135 

0.069 

Relaxation  time 

e  x  10^  sec  ^  5.48 

6.93 

10.29 

13.47 

' 


RELAXATION  TIME  (  sec  ) 


FIGURE  A- 2 


SHEAR  STRESS  AT  WALL 
(  psf  ) 
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Table  A-IV  represents  the  predicted  friction  factors 


for  the  1 " 

and  2.75"  pipes 

using  the  relaxation 

time  estimated  in 

Figure  A- 

2.  The  predi cted 

friction  factors  are  plotted  in 

Fi  gure 

IV- 1 . 

TABLE  A-IV 

PREDICTED  DRAG  REDUCTION 

1" 

Pipe 

2.75" 

Pi  pe 

Reynolds 

fri ction 

friction 

fri cti on 

fri  cti  oi 

number 

factor 

veloci  ty 

factor 

velocity 

(ft/sec) 

(ft/sec 

* 

* 

Nd 

Re 

f 

U 

f 

U 

5,000 

0.00873 

0.0475 

0.00952 

0.0180 

10,000 

0.00628 

0.0807 

0.00748 

0.0320 

20,000 

0.00438 

0,1340 

0.00575 

0.0562 

50,000 

0.00264 

0.2610 

0.00382 

0.1140 

100,000 

0.00166 

0.4150 

0.00272 

0.1930 

200,000 

0.00107 

0,6680 

0.00185 

0.3180 

• 
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iv)  B  Function 

Obtained  from  Velocity  Profiles 

Table 

A-V  represents  the  intercept  values  of  the 

corrected 

semi -1 ogari thmi c 

velocity  profiles  in  Figure  IV-17.  These 

values  of 

B  for  the  0.01% 

solutions  are  plotted  versus  the  friction 

velocity 

in  Figure  IV- 18. 

TABLE  A-V 

VALUES  OF  B  FUNCTION 

Run 

* 

U 

B(e+) 

Number 

friction  velocity 

(ft/sec) 

6 

0  o 1 890 

11.0 

16 

0.1832 

12.5 

10 

0.1830 

14.7 

12b 

Or  3606 

19.2 

12a 

0.3462 

%  21.0 

APPENDIX  B 


B-l 


CALIBRATION 

i )  Calibration  of  the  Magnetic  Flowmeter 

The  bulk  flow  rate  was  measured  by  a  2"  Foxboro,  Model  696, 
magnetic  flowmeter  connected  to  a  Hewlett  Packard  strip  chart  recorder. 
From  a  previous  study  (6),  calibration  curves  for  water  were  obtained 
from  a  weighing  procedure  for  two  sensitivities  in  order  to  cover  a 
range  of  flow  rates  from  0-32  and  0-210  USgal/min.  Figure  B-l 
shows  the  high  and  low  sensitivity  curves  of  flow  rate  versus  voltage 
obtained  from  the  previous  study.  The  data  points  presented  on  the 
figure  serve  as  check  points.  These  data  were  taken  during  the  course 
of  this  work  and  verify  the  calibration  curves. 

i i )  Check  of  Velocity  Profiles 

For  every  photographic  run,  a  check  of  the  velocity  profile 
against  the  flowmeter  reading  was  obtained  by  integrating  the  velocity 
profile.  The  percent  difference  between  the  bulk  velocities  obtained 
from  the  integrated  profiles  and  the  recorded  bulk  velocities  are 
shown  on  Table  B-I.  A  +  2%  average  difference  can  be  found  (except 
for  run  12  which  was  -4 %)  confirming  good  agreement  for  the  measured 
profi les . 
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voltage  (volts  ) 


4 

6 

9 

10 

11 

12 

13 

14 

15 

16 

17 
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TABLE  B-I 


COMPARISON 

OF  RECORDED 

BULK  VELOCITY 

WITH  INTEGRATED  BULK  VELOCITY 

BULK  VELOCITY 
(ft/sec) 

PERCENT 

DIFFERENCE 

Recorded 

Integrated 

4.825 

4.802 

+0.48 

4.428 

4.428 

-1.76 

4.975 

5.084 

-2.15 

4.980 

5.080 

-1.97 

13.26 

13.672 

-2.49 

11.62 

12.128 

-4.18 

22.02 

21.616 

+1.87 

33.40 

33.64 

+0.72 

4.996 

5.142 

+2.79 

5.068 

5.09 

+0.43 

8.77 

8.668 

+1.08 

■ 
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i i i )  Axis  of  the  Digitizer 

As  already  mentioned,  the  full  scale  movement  of  the  digitizer 
cross  hair  for  both  axes  of  the  digitizer  was  10  volts.  In  order  to 
transform  the  digitizer  readings  (in  mvolts)  into  length  scales, 
calibration  functions  for  both  axes  were  obtained,  mainly  in  the  limits 
of  the  projected  frames  by  taking  readings  of  a  10  cm  ruler. 

Figure  B-2  shows  the  voltage  corresponding  to  10  cm  move¬ 
ment  in  both  the  X  and  Y  directions  for  positions  spanning  the  digitizer 
screen.  Calibration  points  in  the  X-direction  have  been  determined  for 
various  Y  positions  and  similarly  Y-direction  readings  have  been  taken 
for  various  X-positions.  The  data  show  the  response  is  constant  over 
the  entire  surface  of  the  screen.  Before  analyzing  a  series  of  photo¬ 
graphs  these  data  were  used  as  a  check  to  ascertain  that  the  digitizer 
electronics  were  stable  and  operating  properly, 

iv)  Projected  Axial  Magnification 

As  mentioned  in  Chapter  III,  several  pictures  of  a  ruler 
divided  into  1/10  cm  increments  were  taken  and  projected  on  the 
digitizer  screen  where  the  projected  length  was  measured.  Because 
difficulties  were  encountered  in  disconnecting  the  optical  box  for 
every  run,  in  order  to  take  a  picture  of  the  ruler  inside  the  pipe, 
the  pictures  were  taken  with  the  ruler  outside  the  pipe.  To  assure 
the  validity  of  this  procedure,  the  projected  axial  magnification 
was  obtained  from  two  different  rulers,  one  inside  the  pipe,  the 
other  outside.  As  shown  in  Table  B-II  for  run  15  as  an  example,  the 
percentage  difference  in  magnification  for  the  photographs  of  rulers 
inside  and  outside  the  pipe  demonstrates  the  correctness  of  the  procedure. 
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FIGURE  B-  2 
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Very  careful  attention  was  given  to  the  focus  of  the  ruler, 
particularly  with  the  1.9  lens  assembly  which  had  a  relatively  large 
depth  of  focus.  Table  B-III  shows  results  for  run  15  where  magnifications 
were  obtained  for  different  focal  distances  (the  ruler  still  in  focus). 

For  a  maximum  deviation  of  30  microns  from  the  apparent  point  of 
focus,  the  magnification  could  be  off  by  1  %.  This  indicates  little 
variability  in  magnification  owing  to  focus  of  the  camera. 


TABLE  B-II 

AXIAL  MAGNIFICATION  DATA  FOR  RUN  15 


Projected  axial  magnification 


i nside 

ruler  outside 

31 .2 

31 .0 

31.7 

31 .5 

31 .3 

31  .5 

31  .5 

31.7 

31.5 

31  .7 

31 .2 

31.5 

31 .7 

31,5 

31 .9 

31 .7 

31 .8 

31 .8 

31.1 

— 

31.49  cm/  cm 

31 . 54  cm/  cm 

Percent  difference  is  0.159% 


TABLE  B-III 

AXIAL  MAGNIFICATION  WITH  CHANGE  OF  FOCUS 

For  the  camera  raised  higher  than  the  focal  point  for  run  15. 

focal  point  10  microns  20  microns  30  microns 

31 .50 
31  .50 
31.50 


31 .30 
31.54 


31.50 

31.30 

31.46 


31.12 
31  .34 
31  .24 


■ 
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v )  Radial  Magnification 

As  in  the  case  of  defining  the  axial  magnification,  projected 
radial  lengths  were  measured  to  determine  the  radial  magnification. 

As  explained  in  Chapter  III,  Section  1 1 1-6,  photographs  of  the  tip  of 
a  needle  were  taken  at  known  radial  positions.  Figures  B-3  and 
B-4  show,  as  an  example  for  run  6,  the  curves  obtained  for  high  and 
low  magnification  of  actual  distance  from  the  wall  versus  the  distance 
measured  on  the  digitizer  screen.  It  can  be  noted  that  the  curve  for 
high  magnification  crosses  the  origin  and  that  pipe  curvature  makes 
the  magnification  vary.  A  third  order  Lagrangian  interpolating  poly¬ 
nomial  was  used  in  the  computer  program  to  represent  these  curves. 

Data  obtained  for  calibrations  for  each  of  these  runs  in 
this  work  have  been  omitted  from  this  thesis  because  they  are  of  no 
significance  once  the  projecting  apparatus  has  been  moved.  The 
axial  and  radial  magnifications  can  be  determined  directly  from  the 
photographs  which  have  been  filed. 
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APPENDIX  C 
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ANALYSIS  OF  ERRORS 


i )  Confidence  Interval  for  Time  Average  Velocities  and  Variances 


The  mean  values  and  variances  of  the  axial  and  radial  velo¬ 
cities  were  obtained  from  measurements  of  instantaneous  velocities. 

Because  the  population  from  which  the  sample  was  drawn  has  unknown  mean 

o 

(y)  and  unknown  variance  (a  ),  the  100(1  -  a)%  confidence  interval  for 
the  mean  value  in  which  y  is  contained,  if  30  or  more  observations  are 
taken,  is 


u  +  Z  (-24 

- 


C-l 


and  the  100(1 


confidence  interval  for  a 


(n-U^  2  (n-l)S* 

- -  1  a  <_  -2 - — 

Xn-1  ,a/2  Xn-1  ,1 -a/2 


C-2 


where  S  is  the  variance  of  a  sample  of  n  observations  from  a  normal 
m  r 

2  * 

distribution  N(y,a  )  . 

As  already  mentioned  and  tabulated  in  Appendix  H  most  of  the 
mean  velocities  and  variances  were  obtained  from  a  large  sample  (larger 
than  80).  In  order  to  show  the  confidence  limit  associated  with  the 
data  of  this  work  an  example  is  done  for  the  case  where  the  number  of 
observations  was  small.  This  gives  the  largest  confidence  interval. 


"Introductory  Engineering  Statistics",  by  Guttman,  I.  and  Wilks,  S.S., 
John  Wi ley ,  N. Y.  ,  1965. 
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For  run  6  at  a  radial  position  near  the  wall  27  observations  were 
recorded  and  the  sample  mean  velocity  and  variance  (in  the  axial 
direction)  obtained  were 

u  =  1 . 958  ft/sec 

X 

S  =  0.568  ft/sec 
m 

Thus  the  95%  confidence  interval  for  the  mean  velocity  is 

+  0.318  S 
—  m 

and  for  the  variance 


0.818  S  <a<  1 .373  S 

m -  m 

Thus  the  time  averaged  velocities  are  generally  estimated  to  within 
less  than  5%  at  95%  confidence.  The  confidence  interval  for  variance 
is  somewhat  larger. 

i i )  Experimental  Errors  in  Reading  the  Streaks 

The  mechanical  errors,  e,  involved  in  reading  the  streaks 
from  the  digitizer  table  can  be  expressed  as 


where  Xm  is  the  measured  length  of  a  streak 
is  the  true  length  of  a  streak. 

If  no  intrinsic  error  is  assumed  and  the  mechanical  errors  are  normally 
distributed,  then  the  net  effect  on  the  measured  mean  value  is  zero  and 
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In  the  case  of  the  sample  variance,  the  mechanical  errors  are  additive. 
The  measured  sample  variance  can  be  expressed  as 

S2  =  ~T  l  (x  '  X) 2  C-5 

m  n-1  ^  '  m 

Substituting  Equation  C-3  in  the  last  equation  and  manipulating  the 
terms  the  following  expression  is  obtained 


S2  =  J- 
m  n-1 


n  ? 

H  + 


n-1 


I  (xt  -  x)i  +  fPT  f  2e-i (xt  -  c'6 


The  first  term  in  this  equation  can  be  estimated  from  repeated  readings 
of  a  known  length  with  digitizer  apparatus.  The  last  term  of  Equation 
C-6  represents  the  covariance  between  the  error  and  (X^  -  X)  and  is 
zero  because  random  mechanical  errors  cannot  be  correlated  with  turbu¬ 
lent  fluctuation  (X^  -  X) . 

Thus  Equation  C-6  can  be  expressed  as 

?  2  ? 

S  =  S  +  C-7 

m  e  T 


or  rearranging  the  equation 


To  estimate  Se,  Figure  C-l  gives  the  variance  obtained  from  reading  a 
given  streak  many  times  with  the  digitizer  as  a  function  of  the  length 
of  the  streak.  These  readings  show  scatter  or  variance  because  of  the 
digitizer  electronics  error  invol ved> because  of  the  blurry  ends  of  the 
streaks,  and  human  error  in  positioning  the  cross  hair.  The  error 
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increases  slightly  with  length  because  the  longer  streaks  have  a  larger 
tail  or  blurried  region  and  it  is  therefore  more  difficult  to  locate  the 
beginning  or  end  of  the  streak. 

As  an  example,  run  14  was  chosen  because  of  the  measured  low 
intensities.  The  calculated  sample  variances  obtained  at  y  =  0.0119  are 

Smx  =  2.76  ft/sec 

S  =  0.239  ft/sec  for  n  =  87 
mr 

as  found  in  Appendix  H. 

At  this  position  the  corresponding  average  actual  lengths 
measured  were 

ax  =  0.100"  and  Ar  *  C.000" 

Taking  the  axial  magnification  (31.4)  and  using  the  radial  calibration 
curve  both  lengths  can  be  transferred  into  projected  lengths  on  the 
digitizer  screen 

ax  =  3.05"  and  Ar  -  0.00" 

Using  Figure  C-l  ,S  can  be  found  as 

S  =  0.0084"  and  =  0.010" 
ex  er 

Or  in  terms  of  velocities,  using  the  time  scales  in  Table  C-I 

S  =  0.0398  ft/sec,  Ser  =  0.1485  ft/sec 
ex 


-  S 

tx  mx 


S.  =  0.785  S 
tr  mr 


Plugging  into  Equation  C-8 
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This  calculation  shows  that  errors  in  the  calculated  axial  intensity 
owing  to  measuring  errors  are  entirely  negligible.  For  the  radial 
intensity  however,  the  calculated  value  can  be  as  much  as  21.5%  higher 
than  the  true  value  owing  to  the  mechanical  error.  This  is  the  extreme 
case  and  generally  for  other  runs  or  other  radial  positions  the  percent 
error  will  be  much  smaller. 

i i i )  Location  of  Wall  and  Tracers 

From  Figure  1 1 1 -6b  it  can  be  observed  that  the  wall  location 
is  masked  by  a  shadow  region.  Along  with  uncertainty  because  of 
wall  curvature  it  was  estimated  from  a  projected  photograph  that  the 
wall  could  be  located  to  within  +  0.0010".  The  smallest  reading 
recorded  for  the  radial  position  near  the  wall  is  for  run  6  and  is 
0.0040".  This  radial  position  can  then  be  25%  in  error,  but  represents 
the  worst  case  encountered.  Furthermore  the  error  of  +_  0,001"  in  the 
locating  the  wall  is  of  the  same  order  as  the  length  corresponding 
to  non-zero  average  radial  velocities  reported  in  Appendix  H. 

i v)  Error  in  Intensity 

a )  air  bubbles 

As  already  discussed  in  Seyer's  thesis  (41  ) ,  the  influence  of 
the  air  bubbles  on  the  intensity  measurements  will  be  negligible  if 
their  size  is  smaller  than  the  scale  of  the  energy  containing  eddies 
(approximately  0.006").  Because  in  this  work  the  average  diameter  of 
the  bubbles  was  found  to  be  0.002",  no  influence  is  expected. 

b )  uncertainty  in  radial  position 


The  fact  that  velocities  were  taken  over  a  band  of  positions 


' 
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rather  than  at  a  single  radial  position  causes  the  calculated  intensities 
to  be  too  high.  This  occurs  because  of  the  mean  velocity  gradient  and 
so  velocities  of  streaks  at  the  bottom  of  the  band  (nearest  to  wall)  are 
on  the  average  slightly  smaller  than  velocities  of  streaks  near  the  top. 
Thus  there  is  a  fluctuation  in  velocity  owing  to  variations  in  the 
radial  position  of  the  streaks.  The  magnitude  of  the  error  depends  on 
the  velocity  gradient  as  well  as  the  width  of  the  band. 

In  reading  the  photographs  the  radial  positions  of  interest 
were  determined  by  lines  drawn  on  a  clear  plastic  sheet  which  was 
placed  on  the  digitizer  screen.  The  situation  is  illustrated  by 
Figure  1 1 1-7  where  the  limits  of  the  band  correspond  to  lines  on  the 
plastic  sheet.  For  positions  near  the  wall  the  spacing  of  the  lines 
on  the  sheet  was  approximately  .25"  and  was  about  1"  for  positions 
near  the  centre.  Fore  the  high  magnification  portion  of  run  6  the 
.25"  band  on  the  plastic  sheet  corresponds  to  a  band  ,0024  inches 
wide  in  the  one  inch  pipe, 

A  conservative  estimate  of  the  error  in  the  calculated  inten¬ 
sity  can  be  obtained  by  assuming  the  streaks  are  located  either  at  the 
center  or  top  or  bottom  of  the  band.  In  fact  there  will  be  a  distri¬ 
bution  of  positions  which  peaks  sharply  at  the  center  of  the  band.  In 
the  following  sketch  the  positions  y  +  Ay  and  y  -  Ay  represent  the 
limits  of  the  band  for  a  radial  position  y.  Suppose  that 

- - - — - - —  y  +  Ay,  u  +  Ail 

acceptable  streaks  J - - - - -  y,  u 

k - _____  y  -  Ay,  u  -  alT 

n  streaks  are  observed  in  the  band.  Roughly  n/2  streaks  will  have  a 
velocity  characteristic  of  position  y,  while  n/4  streaks  will  have 
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velocities  reflecting  the  top  position  and  n/4  reflecting  bottom, 
Thus  the  variance  owing  to  velocity  gradient  will  be 


S2  =  ’ 

u  n-1 


I(U  -  U)f 


n/2 


-  —\2 


3/4n 


n-1 


I  (U  -  IT)  +  l  (u  -  (u  +  au")  ) f 
1=1  n/2  1 


+  l  (U  +  (U  -  AU  )  )  ^ 

3/4n  1 


1  (  a — a  ^ 

*  2"  (Au) 
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If  Sm  and  are  the  measured  and  true  variances  then 


2  2  2 

S  =  87  + 
m  t  u 


or 


S  Z 

(A4  c-io 

m 


This  equation  shows  the  influence  of  the  velocity  gradient  on  the  true 
i ntensi ty . 

As  an  example  we  will  consider  the  first  3  data  points  in 
Table  H- 1-3  for  Run  6.  Consideration  of  the  velocity  profiles  shows 
this  case  will  have  the  largest  error  owing  to  the  above  considerations. 
In  the  following  table  the  appropriate  data  for  use  of  equation  C-10  is 


tabulated.  The  Au 
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No. 

y 

i  nch 

AU 

ft/sec 

S 

m 

ft/sec 

S./S 
t  m 

1 

.0040 

.18 

.567 

.95 

2 

.0073 

.13 

.508 

.97 

3 

.0102 

.04 

.535 

.97 

for  a  band  width  of  .0024  inches  have  been  determined  directly  from 
the  slopes  of  a  plot  of  the  velocity  profile  at  the  indicated  radial 
positions.  The  measured  intensity  is  obtained  directly  from  Table 
H- 1-3  by  multiplying  I  by  the  maximum  velocity. 

X 

It  is  evident  from  the  tabulated  values  of  S./S  there  is 

t  m 

little  difference  between  the  measured  and  true  intensities.  For  the 
radial  position  nearest  to  the  wall,  where  the  velocity  gradient  is 
largest,  the  error  would  be  approximately  5%.  In  view  of  the  conser¬ 
vative  nature  of  this  calculation  it  is  concluded  that  for  all  runs  in 
this  work  error  in  intensity  because  of  velocity  gradient  is  negligible. 

v)  Time  Scale  Relative  to  the  Dissipative  Frequencies 

Table  C-I  compares  the  time  scales  used  for  each  photographic 
run  and  the  estimated  frequency  of  maximum  energy  dissipation  for 
organic  solvents  as  found  by  Patterson  and  Zakin  (32)  for  1  and  2" 
pipes.  It  can  be  observed  that  for  all  the  runs,  except  run  14,  the 
time  scale  used  are  lower  than  the  dissipative  frequencies.  A  true 
indication  of  the  turbulence  intensity,  (containing  the  contribution 
from  all  significant  frequencies)  was  therefore  obtained. 
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TABLE  C-l 

TIME  SCALE  VERSUS  DISSIPATIVE  FREQUENCIES 


RUN  NUMBER 

BULK  VELOCITY 

N 

TIME  SCALE 

ft/sec 

max, 

sec" 

sec  ^ 

4 

4.8 

250 

596.6 

5 

3.1 

130 

596.6 

6 

4.4 

220 

596.6 

7 

5.0 

260 

596.6 

9 

5.0 

260 

1194 

10 

5.0 

260 

1197 

11 

13.2 

1000 

1787.5 

12 

11.6 

800 

1787.5 

13 

22.0 

1900 

1787.5 

14 

33.4 

~4000 

1787.5 

15 

5.0 

260 

895 

16 

5.0 

260 

895 

17 

8.8 

600 

895 

. 
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SAMPLE  CALCULATIONS 

i )  Velocities 

The  sample  calculation  will  be  made  for  water  run  11  at  a 
radial  position  y  =  0.0906  inch  from  data  tabulated  in  Appendix  H. 


Conditions  of  Operation 


Temperature 

Flow 

Manometer  Reading 
Length  between  taps 
Diameter  of  tube 
Rotation  of  disk 
Number  of  slits 


70°F 

32.4  gal/min 
6.1  in  of  Hg 
11.0  ft 


1.0  in 


3575  rpm 


60  slit-spoke/rev 


Readings  from  the  Digitizer  of  Streak  Relative  to  Pipe  Wall 

Referring  to  Figure  1 1 1-7  the  following  points  correspond 
to  the  reading  of  the  streak  and  spoke  and  the  wall.  The  points  are 
in  volts  and  relative  to  the  digitizer  axis  as  zero. 


Y-|  =  6.287  volts 
Y2  =  6.221  volts 
Y^  =  8.077  volts 
Y^  =  8.052  volts 
With  calibration  constants 

k  =  0,1718  volt/cm 
x 


X  -  5,401  volts 
=  3.363  volts 
=  7.552  volts 
X^  =  4.972  volts 
as  shown  in  Figure  B-2 

k  =  0.1811  volt/cm 


these  readings  can  be  transformed  into  lengths  on  the  projection  screen. 
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Yn  =  13.668  inches 


X,  =  12.377  inches 


=  13.524  inches 
=  17.559  inches 
Y^  =  17.505  inches 


=  7.707  inches 

=  17.306  inches 
X^  =  11 .394  i nches 


With  these  numbers,  from  the  geometry  of  the  figure  the  axial  component 
of  the  streak,  ax,  parallel  to  the  wall  is  found  to  be 


Ax  =  4.72  inches 


The  measured  axial  magnification  is  53.18  and  therefore  the  true 
length  of  the  streak  is  4.72/53.18  =  .0887  inch. 

From  the  rotational  speed  of  the  chopping  wheel  and  number  of 
slit-spoke  of  the  wheel  the  time  scale  can  be  found 

T  =  (60  sec/min)/(30  division/rev)(3575  rpm)  =  0.00056  sec 
and  the  instantaneous  velocity  is 


u  =  ax/T (53.18)  =  13.184  ft/sec 

X 


Similarly,  the  radial  instantaneous  velocity  was  found  to  be  (using  the 
radial  calibration  curve,  as  shown  for  run  6  as  an  example,  in  Figures 
B-III,  B- IV ) 


u  =  -0.330  ft/sec 
r 


From  the  computer  the  average  of  all  these  readings  gives  for  the 
mean  velocities 


u  =  13.400  ft/sec  and  u  =  +0.219  ft/sec 

X  1 


and  the  variances  S  =  1.022  ft/sec  and  S  -  0.47  ft/sec. 

X  1 


' 


. 
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i i )  Flow  Calculations  for  Water 

The  pressure  drop  between  two  pressure  taps  can  be  calculated 
from  the  manometer  reading 

AP  ■  I.6.,l,|,(62.3j.jl2.55)  .  397,5  psf 

The  Reynolds  number 

NRe  =  <Irx>D/v 

NRe  =  Q  D  (8.337 )/ (pAv)  =  104299 

The  shear  stress  at  the  wall 

t. ,  =  DAP/4L  =  0.753  psf 
w 

The  friction  factor 

f  =  2tii/(p<u  >2)  =  0.00442 

W  X 

The  friction  velocity 

U*  =  A  g  /p  =  0.624  ft/sec 
w  c 

The  dimensionless  distance  from  the  wall 

Y+  =  y  U*/v  =  444.5 

The  dimensionless  velocity 

U+  =  u  /U*  =  2.148 

X 


. 
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iii)  Calculations  of  K1  ,  n'  and 

A  sample  calculation  is  made  for  run  13  as  an  example  for 
Separan  AP-30  0.10%  and  these  data  can  be  found  in  Appendix  H, 


Flow  54.00  USgal/min 

Manometer  reading  7.95  inches  of  Hg 

Length  between  taps  22.00  ft 

From  these 

AP  =  519.38  psf 

t  =  0.491  psf 
w 

The  constants  of  viscometric  equation  are  (see  Appendix  A) 
A  =  -6.540  B  =  0.5989  C  =  0.00599 


so  that 


In 


8<u  >/D 
x 


-0.5989  +  / (0.5989)2  -  4 ( 0 . 0060 ) ( - 6 . 540  -  In  0.491) 

2(0 .00599 ) 


In  3<u  >/D  =  8.95 

A 

8<u>/D  =  7.71  x  103  sec"1 

X 

n'  =  0.5989  +  2(0.00599)(8.95)  =  0.705 
K'  =  0.491/(7.71  x  103)0,705  =  8.93  x  10 
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i v)  Determination  of  e  from  Friction  Data 

The  relaxation  time  found  from  friction  data  in  the  1"  pipe 
is  tabulated  in  Table  A-3.  From  Appendix  I  the  friction  factor  is 


f  =  0.00177 
and  NRe  =  93369 

Putting  the  constants  in  Equation  1 1 1 - 1 6  the  following  expression  is 
obtained 

—  =  1.74(1  -  sp2  In  (NRe/T)  +  0.707(1  -  51)2[B(e+) 

vf 

-  2.56)]  -  — 

Sf 


From  the  data 

/f  =  0.0421 


l//f  =  23.753 


NRg/f  =  3931  ln(NRe/f)  =  8.2764 

For  the  first  trial  assuming  e,- j  =  0,  we  have,  solving  the  above  equation 
for  B, 

B(e+)  =  (23.753  -  14,401  +  1.810  +  2.121)/0.707  =  18.79 


Thus  using  Equation  11-12,  the  first  estimate  of  c-j  is 

Y|  -2.46  In  Y|  =  18.79 
y|  =  26.9  and  5-j  =  0.0194 


and 


I 
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For  the  second  trial  ^  =  0.0194  is  used  as  an  approximation  and 

B(e+)  =  20.25 
Y|  =  28 

^  =  0.0202 

The  final  B(e  )  which  satisfies  the  above  equation  is 
B(e+)  =  20.30.  From  the  Seyer's  B  function  (see  Figure  1 1 -3 ) 
e+  =  9.25  and  from  the  experimental  values  of 

U*  =  0.4002  ft/sec 

x  =  0.3096  psf 
w 

6  =  6+  v/U  2  =  6,93  x  10~4  sec 
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TURBULENT  SHEAR  STRESS  AND  KINETIC  ENERGY  DATA 

In  order  to  provide  additional  support  for  the  turbulence 
measurements,  it  is  possible  to  estimate  the  covariance  of  the 
instantaneous  velocities  as  well  as  the  axial  kinetic  energy  of  the 
turbulence 

1  n 

S  =  Y  y  [(u  *  u  )(u  -  u  )].  E-l 

xr  n-1  ^LVx  x/vr  r;ji 

which  is  proportional  to  the  Reynolds  shearing  stress.  Similarly  for 
the  kinetic  energy  in  axial  direction 

S  =  '  V  [(u  -  u  )dl.  E-2 

xx  n-1  |  LV  x  x'  Ji 

Sample  calculations  were  done  for  two  radial  positions  for 
water  run  11  in  the  1"  pipe  because  of  the  large  number  of  observations 
available.  From  Appendix  Kin  Book  II  the  instantaneous  velocities 
can  be  found  for  positions  y  -  0.462  inch  and  y  =  0.387,  and  the 
mean  velocities  can  be  found  in  Appendix  H. 

Table  E-I  represents  the  calculated  Reynold  stress  a  id 
the  kinetic  energy.  These  terms  are  plotted  on  Figure  E-l  and 
compared  to  Laufer's  data. 
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TABLE  E-I 


TURBULENCE 

PROPERTIES 

Radial  position  -y  (inch) 

0.462 

0.387 

Number  of  observations 

210 

192 

Vr 

-0.0483 

0.0098 

— 1 - 1 - 

u  u 

X  X 

0.6859 

0.7075 

'  '  *p 

Vr/U 

-0.1242 

0.0252 

1  1  *o 

Vx/U 

1.763 

1.819 

Three  of  the  four  points  are  in  excellent  agreement  with  the  expected 
values.  The  low  point  for  shear  stress  is  in  fact  in  reasonable 
agreement  in  view  of  the  extreme  difficulty  in  estimating  a  covariance. 
Reference  to  statistical  tables  shows  that  with  these  points  the  95% 
confidence  interval  about  the  lower  point  includes  the  known  curve 


of  shear  stress. 


. 
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TUBBULENT  KINETIC  ENERGY 
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FIGURE  E-1 


REYNOLDS  STRESS  AND  KINETIC  ENERGY 


TURBULENT  SHEAR  STRESS 
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Rq  correlation  factor 

The  R^  correlation  factor  given  in  Table  F-I  and  defined  in 

Equation  II-2was  calculated  for  all  the  radial  positions  of  all  the 

runs  in  this  work  and  also  for  some  of  Seyer's  data  (4l)»  RQ  is 

simply  the  ratio  of  number  of  positive  radial  velocities  to  the  number 

of  negative  radial  velocities.  To  calculate  R  raw  data,  that  is  each 

o 

instantaneous  radial  velocity  in  Appendix  L,  has  been  adjusted 
slightly  by  requiring  that  the  mean  radial  velocity  be  exactly  zero. 

A1 ternati vely  the  non-zero  mean  radial  velocity  has  been  subtracted 
from  each  instantaneous  velocity.  This  effectively  removes  bias  in 
the  data  owing  to  small  error  in  locating  the  wall.  This  has  the 
effect  of  changing  the  relative  number  of  negative  and  positive 
fluctuations;  but  only  insignificantly. 


. 
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TABLE  F-I 


Run 


n 

6=y/R 

R 

0 

No .  n 

C=y/R 

R 

30 

0.0302 

1.14 

6  27 

0.0079 

1. 

74 

0.0476 

0.92 

39 

0.0146 

1. 

63 

0.0617 

1 .36 

33 

0.0178 

2. 

61 

0.0790 

1 .48 

53 

0.0204 

2. 

77 

0.0949 

1  .43 

64 

0.0311 

3. 

79 

0.1092 

1  .30 

51 

0.0389 

1. 

75 

0.1250 

1.50 

79 

0.0467 

2. 

72 

0.1450 

1  .49 

60 

0.0540 

1. 

67 

0.1575 

1.16 

67 

0.0614 

1  . 

70 

0.1735 

1.00 

67 

0.0689 

1. 

62 

0.1890 

0.83 

47 

0.0743 

1. 

47 

0.2050 

1  .86 

32 

0.0783 

1  . 

50 

0.2205 

1  .38 

29 

0.0873 

2. 

63 

0.2390 

1.10 

78 

0.0966 

1 . 

15 

0.0141 

0.67 

79 

0.1756 

1. 

22 

0.0203 

1.75 

37 

0.02543 

1. 

26 

0.0316 

1.60 

34 

0.333 

1  . 

37 

0.0375 

0.54 

23 

0.412 

2. 

37 

0.0473 

0.95 

32 

0.0559 

0.883 

11 

0.0638 

o 

25 

79 

67 

12 

27 

83 

16 

40 

79 

23 

24 

29 

22 

23 

19 

64 

43 

29 


• 

' 

F-3 


(Table 

F-I  - 

Continued) 

Run 

Run 

No . 

n 

s=y/R 

R 

0 

Mo .  n 

e=y/R 

R 

0 

9 

22 

0.012 

2.66 

10  43 

0.016 

2.91 

31 

0.021 

1.58 

62 

0.022 

1.14 

34 

0.029 

0.70 

48 

0.032 

3.00 

49 

0.034 

0.75 

58 

C.040 

2.62 

68 

0.044 

0.79 

79 

0.048 

1.39 

37 

0.052 

0.76 

79 

0.062 

1  .92 

44 

0.059 

1.44 

79 

0.078 

1.83 

51 

0.067 

2.00 

79 

0.094 

2.16 

52 

0.075 

1.26 

79 

0.124 

2.16 

43 

0.084 

1.05 

79 

0.050 

1  .92 

59 

0.093 

1.36 

79 

0.086 

1.25 

75 

0.132 

1  .42 

79 

0.130 

1.34 

73 

0.111 

1.43 

79 

0.206 

1.08 

73 

0.149 

0.59 

79 

0c  286 

1.55 

51 

0.169 

1.22 

79 

0.366 

1 .82 

27 

0.034 

0.50 

79 

0.460 

1  .72 

55 

0.058 

0.83 

79 

0.534 

1  .82 

83 

0.073 

1.13 

55 

0.602 

0.72 

97 

0.112 

0.94 

55 

0.682 

1.03 

95 

0.190 

0.86 

52 

0.832 

2.25 

70 

0.308 

0.71 

56 

0.464 

1  .00 

39 

0.621 

0.44 

21 

0.855 

1  .63 

• 

• 

F-4 


(Table 

F-I  - 

Conti nued) 

Run 

No. 

n 

C=y/R 

Ro 

Run 

No .  n 

5=y/R 

Ro 

11 

18 

0.032 

0.64 

12  15 

0.019 

2.00 

20 

0.049 

1  .86 

18 

0.028 

3.50 

64 

0.064 

1.21 

28 

0.037 

4.60 

79 

0.104 

1.14 

55 

0.044 

3.23 

69 

0.145 

1.03 

90 

0.040 

2.85 

79 

0.181 

1  .85 

107 

0.078 

2.34 

79 

0.224 

2.00 

13 

0.046 

1.17 

79 

0.384 

1.17 

28 

0.062 

1.00 

79 

0.460 

1  .26 

76 

0.094 

1.11 

79 

0.614 

0.49 

100 

0.130 

1.18 

79 

0.774 

1  .08 

100 

0.212 

1.56 

79 

0.924 

2.04 

107 

0.366 

1.68 

94 

0.181 

1  .19 

97 

0.524 

0.90 

123 

0.698 

2.32 

75 

0.684 

0.97 

110 

0.843 

1  .62 

85 

0.842 

1.18 

99 

0.384 

1.25 

193 

0.774 

0.77 

221 

0.924 

1.11 

. 

, 

F-5 


(Table 

F-I  - 

Conti nued) 

Run 

No. 

n 

?=y/R 

13 

72 

0.019 

57 

0.028 

63 

0.036 

64 

0.052 

63 

0.666 

59 

0.082 

16 

0.012 

46 

0.048 

103 

0.090 

103 

0.130 

95 

0.212 

95 

0.290 

95 

0.368 

87 

0.600 

87 

0.760 

R 

0 

Run 

Mo. 

3.00 

14 

5.35 

6.88 

8.14 

11.60 

10.80 

4.33 

1  .88 

4.15 

5.06 

4.00 

6.31 

7.64 

0.67 

16.40 

n 

?=y/R 

Ro 

99 

0.024 

7.25 

102 

0.064 

6.29 

98 

0.098 

3.66 

103 

0.142 

3.90 

106 

0.258 

9.60 

102 

0.422 

9.20 

112 

0.498 

3.00 

100 

0.576 

7.33 

98 

0.734 

11.25 

100 

0.948 

0.22 

. 

. 


F-6 


(Table  F-I  -  Continued) 


n 

c=y/R 

R 

0 

Run 

Mo .  n 

5=y/R 

R 

0 

11 

0.008 

0.83 

16  89 

0.018 

1.87 

59 

0.022 

0.69 

95 

0.039 

1 .07 

78 

0.040 

0.90 

95 

0.075 

1 .57 

95 

0.056 

0.90 

95 

0.146 

1.38 

95 

0.069 

0.98 

95 

0.218 

1.44 

95 

0.111 

0.94 

95 

0.367 

1.57 

95 

0.149 

0.98 

95 

0.511 

1  .97 

95 

0.226 

1.16 

75 

0.657 

0.66 

95 

0.298 

1  .32 

70 

0.808 

2.04 

45 

0.441 

0.96 

90 

0.945 

2.75 

95 

0.510 

1.02 

17  22 

0.015 

1 .75 

42 

0.583 

0.83 

60 

0.037 

2.53 

58 

0.728 

0.93 

91 

0.074 

1  ,94 

26 

0.873 

0.63 

95 

0.114 

4.59 

40 

1  .000 

1.00 

87 

0.148 

2.95 

95 

0.219 

2.52 

79 

0.291 

3.39 

84 

0.436 

2.82 

30 

0.580 

0.58 

33 

0.728 

2.00 

35 

0.871 

2.50 

35 

0.944 

2.50 

. 

F-7 


(Table  F-I  -  Continued) 

DATA  FROM  SEYER 1 S  WORK  (41  ) 

Run  Run 


n 

s=y/R 

R 

0 

No . 

n 

5=y/R 

R 

0 

32 

1.67 

4 

59 

0.074 

1.03 

35 

1.00 

59 

0.146 

1.35 

49 

1 .72 

68 

0.290 

0.79 

46 

1 .42 

76 

0.612 

1.11 

61 

0.702 

0.90 

47 

0.929 

1.35 

29 

1.014 

0.45 

5 

113 

0.073 

1.69 

36 

1.25 

57 

0.658 

1  .28 

55 

1.61 

6 

65 

0.084 

0.76 

56 

1.55 

64 

0,154 

1.20 

77 

0.682 

1  .08 

71 

0.300 

0.97 

50 

0.988 

0.82 

69 

0.618 

1.38 

93 

1.06 

46 

0.938 

0.70 

99 

1.15 

7 

65 

0.083 

1.32 

98 

1.65 

64 

0.155 

1  .20 

58 

0.606 

0.66 

71 

0.300 

1  .03 

40 

0.941 

1.22 

69 

0.618 

1  .38 

46 

0.938 

0.71 

APPENDIX  G 


G-l 


MIGRATION  OF  AIR  BUBBLES 

As  already  indicated  in  Chapter  III,  Section  III-l  ,  the 
tiny  air  bubbles  (average  diameter  of  0.002")  used  as  tracers  were 
uniformly  distributed  in  the  calming  section.  On  the  other  hand 
it  was  noticed  (as  suggested  crudely  by  Figure  1 1 X - 5a  9  for  example) 
that  most  bubbles  were  at  a  radial  position  of  c  -  0.11.  Figure 
G-l  presents  for  photographic  runs  (4,9,16,12,13)  the  number  of 
observed  air  bubbles  as  a  function  of  their  radial  position.  It 
can  be  noticed  that  for  most  of  the  runs  a  peak  is  located  at  the 
radial  position  already  mentioned.  For  water  runs  (4.9)  the  peak  is 
sharp  and  represent  approximately  25%  of  all  bubbles  observed  during 
the  photographic  run.  For  the  case  of  the  photographic  runs  with 
Separan  AP-30  solutions,  the  air  bubble  distribution  is  more  uniform 
in  the  core  region  of  the  pipe. 

Even  though  there  is  an  apparent  migration  of  bubbles  to 
the  radial  position  of  0.1,  the  concentration  in  volume  %  is  still 
extremely  small  and  no  influence  on  the  turbulence  is  expected.  The 
curves  show  one  of  the  limitations  of  the  technique  for  velocities 
near  the  wall  is  the  small  number  of  bubbles  that  occur  there  and 
therefore  a  large  number  of  photographs  is  needed  to  get  a  reasonable 
number  of  instantaneous  velocities. 


DISTRIBUTION  OF  AIR  BUBBLES  IN  PIPE 


SNOIlVAdBSaO 


RADIAL  POSITION  (  y/R  ) 


APPENDIX  H 


H-l 


SUMMARY  OF  ALL  RUNS 

In  this  appendix  all  of  the  data  have  been  tabulated  in 
Table  H-I  in  a  format  to  summarize  the  velocity  and  intensity  profiles 
at  all  the  radial  positions  for  each  photographic  run.  The  primary 
data  concerning  the  velocity  measurements  can  be  found  in  Appendix  L 
(in  Book  II)  and  the  flow  conditions  are  tabulated  in  Table  IV-I. 

The  axial  and  radial  turbulent  intensities  tabulated  are  made  dimension¬ 
less  to  the  maximum  velocity. 
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T 

N ' 

NRe 

f 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 


H-2 


TABLE  H-I-l 

RUN  4 
water 


= 

1 1 .8  USgal/min 

w 

ii 

o 

Ol> 

o 

psf 

= 

77.5°F 

* 

U 

=  0.259 

ft/sec 

= 

38123 

‘V 

=  4.80 

ft/sec 

= 

0.0057 

<  u  > 

x  max 

=  5.80 

ft/sec 

n 

y 

Ux 

u 

r 

Y+ 

U+ 

[x 

I 

r 

inch 

ft/sec 

ft/sec 

30 

0.0148 

3.344 

+0.011 

30.29 

12.93 

0.130 

0.028 

74 

0.0239 

3.790 

-0.021 

48.92 

14.65 

0.072 

0.034 

63 

0.0309 

4.021 

+0.037 

63.20 

15.54 

0.104 

0.037 

61 

0.0385 

4.225 

+0.027 

78.78 

16.33 

0.085 

0.045 

77 

0.0467 

4.312 

+0.045 

95.69 

16.64 

0.068 

0.046 

79 

0.0540 

4.376 

+0.034 

00 

o 

o 

16.91 

0.080 

0.046 

75 

0.0620 

4.463 

+0.066 

127.14 

17.25 

0.074 

0.045 

72 

0.0720 

4.735 

+0.083 

148.04 

18.30 

0.072 

0.054 

67 

0.0786 

4.791 

+0.063 

161  .18 

18.52 

0.081 

0.053 

70 

0.0865 

4.591 

+0.021 

177.19 

17.74 

0.089 

0.045 

62 

0.0945 

5.134 

+0.056 

193.63 

19.84 

0.084 

0.053 

47 

0.1025 

4.990 

+0.057 

209.89 

19.29 

0.093 

0.054 

50 

0.1096 

4.809 

+0.004 

224.42 

18.59 

0.071 

0.036 

63 

0.1183 

5.036 

-0.012 

242.37 

19.46 

0.089 

0.043 

H-3 


TABLE  H- 1-2 

RUN  5 
water 


Q 

— 

7.50  USgal/min 

Tw 

=  0.0592 

psf 

T 

= 

69 .0°F 

* 

U 

=  0.1749 

ft/sec 

NRe 

= 

23785 

<V 

-  3.070 

ft/sec 

f 

= 

0.00642 

"ux>max 

-  3.800 

ft/sec 

n 

y 

ux 

Ur 

Y+ 

U+ 

!x 

V 

inch 

ft/sec 

ft/sec 

1 

15 

0.0071 

1  .595 

-0.007 

9.59 

9.12 

0.099 

0.022 

2 

22 

0.0103 

1 .805 

+0.026 

13.92 

10.32 

0.127 

0.019 

3 

26 

0.0158 

2.236 

+0.041 

21  .44 

12.78 

0.108 

0.029 

4 

37 

0.0194 

2.361 

-0.035 

26.27 

13.49 

0.098 

0.033 

5 

37 

0.0236 

2.445 

+0.009 

32.07 

13.98 

0.085 

0.043 

6 

32 

0.0280 

2.559 

+0.006 

37.96 

14.63 

0.067 

0.052 

11  0.0319  2.565  -0.012  43.29  14.66  0.103  0.035 


7 


■ 

Q 

T 

N' 

Re 

f 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


H-4 


TABLE  H- 1-3 


RUN  6 


Sol ution 

of  SEPARAN  AP-30 

at  a  cone, 

.  of  0. 

01% 

= 

10.65  USgal/min 

Tw 

0.0690 

i  psf 

= 

99 .0°F 

* 

U 

0.189 

ft/sec 

= 

35931 

<V  “ 

4.38 

ft/sec 

= 

0.0038 

<ux^max  “ 

5.35 

ft/sec 

n 

y 

inch 

ux 

ft/sec 

u 

r 

ft/sec 

Y+ 

U+ 

I 

X 

I 

r 

27 

0.0040 

1  .958 

0.024 

6.93 

10.35 

0.106 

0.013 

39 

0.0073 

2.540 

0.033 

12.8 

13.42 

0.095 

0.016 

53 

0.0102 

3.022 

0.048 

17.9 

15.97 

0.100 

0.016 

64 

0.0156 

3.290 

0.065 

27.4 

17.39 

0.104 

0.015 

51 

0.0194 

3.527 

0.028 

34.1 

18.62 

0.092 

0.016 

79 

0.0233 

3.609 

0.034 

40.9 

19.05 

0.083 

0.014 

60 

0.0270 

3.878 

0.026 

47.3 

20.45 

0.075 

0.019 

67 

0.0307 

3.888 

0.049 

53.8 

20.52 

0.070 

0.020 

67 

0.0345 

3.962 

0.027 

60.5 

20.93 

0.067 

0.019 

47 

0.0372 

4.101 

0.018 

65.3 

21  .64 

0.053 

0.020 

32 

0.0392 

4.166 

0.061 

68.8 

22.00 

0.088 

0.026 

29 

0.0437 

4.162 

0.023 

76.7 

22.00 

0.091 

0.022 

33 

0.0089 

2.310 

0.062 

15.6 

12.21 

0.122 

0.020 

78 

0.0483 

4.112 

0.001 

84.7 

21  .70 

0.072 

0.019 

79 

0.0878 

4.290 

0.023 

154 

22.62 

0.072 

0.032 

37 

0.1272 

4.538 

0.093 

223 

23.92 

0.067 

0.028 

34 

0.1666 

4.834 

0.019 

292 

25.50 

0.070 

0.036 

23 

0.2060 

4.695 

0.056 

362 

24.78 

0.062 

0.028 

7 

0.2450 

5.110 

0.043 

0.025 

0.024 

8 

0.2850 

5.070 

0.109 

0.052 

0.027 

. 

H-5 


TABLE  H- 1-4 
RUN  7 

Solution  of  SEPARAN  AP-30  at  a  cone,  of  0.10% 


Q 

=  12.1 

USgal/min 

Tw 

=0.075  psf 

T 

=  87°F 

★ 

U 

=0.197 

ft/sec 

N' 

Re 

=  5192 

<5X> 

=4.98  ft/sec 

f 

=  0.00320 

^x'  max 

=  7.74 

ft/sec 

n 

y 

inch 

Gx 

ft/sec 

u 

r 

ft/sec 

Y+ 

U+ 

*x 

I 

r 

1 

40 

0.1648 

5.553 

+0.025 

300 

28.15 

0.235 

0.024 

2 

40 

0.2435 

6.506 

+0.019 

447 

32.99 

0.188 

0.032 

3 

40 

0.3222 

6.940 

0.000 

588 

35.19 

0.195 

0 .055 

4 

40 

0.3616 

7.361 

-0.001 

659 

37.33 

0.291 

0.067 

5 

40 

0.4404 

7.650 

0.025 

816 

38.80 

0.220 

0.046 

. 


D 

H-6 


TABLE  H- 1-5 
RUN  9 


water 


Q 

- 

12,15  USgal/min 

Tw 

=  0.1319  psf 

T 

76. 0°F 

* 

U 

=  0.2612  ft/sec 

N 1 

IN  Re 

= 

42229 

<V 

=  4.975 

ft/sec 

f 

= 

0,0055 

<u  > 

x  max 

=  6.227 

ft/sec 

n 

y 

u 

X 

u 

r 

Y+ 

U+ 

!x 

Tr 

inch 

ft/sec 

ft/sec 

1 

22 

0.0061 

2.889 

+0.037 

13.48 

11.06 

0.137 

0.014 

2 

31 

0.0105 

3.642 

+0.051 

23.30 

13.94 

0,096 

0.024 

3 

34 

0.0145 

3.598 

-0.016 

32.20 

13.78 

0.071 

0.026 

4 

49 

0.0174 

3.741 

-0,036 

38.63 

14.32 

0.075 

0.027 

5 

68 

0.0219 

3.946 

-0.027 

48.51 

15.11 

0,094 

0.034 

6 

37 

0.0259 

3.962 

-0.002 

57.52 

15.17 

0.069 

0.025 

7 

44 

0.0294 

4.111 

+0  o 030 

65.21 

15,72 

0,051 

0.029 

8 

51 

0.0334 

4.055 

+0.077 

74.06 

15.52 

0.056 

0.029 

9 

52 

0.0376 

4.242 

+0.007 

83.35 

16.24 

0.064 

0.028 

10 

43 

0.0418 

4.288 

+0.009 

92.62 

16,42 

0.077 

0.028 

11 

59 

0.0462 

4.426 

+0.024 

102,44 

16.95 

0.076 

0.033 

12 

75 

0.0658 

4.747 

+0.025 

145.88 

18.17 

0.065 

0.029 

13 

73 

0.0551 

4,539 

+0.072 

122.24 

17.38 

0,061 

0.030 

14 

73 

0.0745 

4,732 

-0.042 

165.08 

18.12 

0.061 

0.026 

15 

51 

0.0846 

4.910 

+0.023 

187,63 

18,80 

0.068 

0.033 

16 

27 

0.0172 

3.492 

-0.118 

38,14 

13.37 

0,089 

0.041 

17 

55 

0.0289 

3,838 

+0.024 

64,08 

14.70 

0.085 

0.042 

18 

83 

0.0367 

4.269 

+0,001 

81.37 

16.35 

0.080 

0,043 

19 

97 

0.0562 

4.600 

-0,017 

124,61 

17.61 

0.066 

0.049 

20 

95 

0.0952 

4.796 

+0.004 

211.09 

18,36 

0.066 

0.044 

21 

70 

0.1537 

5.204 

-0,025 

340,80 

19,92 

0.070 

0,039 

22 

56 

0.2317 

5.688 

+0.014 

513.75 

21.78 

0.051 

0,034 

23 

39 

0.3097 

6.035 

-0,076 

686.70 

23.11 

0.051 

0.035 

24 

21 

0.4267 

6.222 

+0,092 

946.12 

23.82 

0.031 

0.026 

*  25 

50 

0.0652 

4.645 

+0.028 

144.56 

17.78 

0,0561 

0.036 

*  taken  separately 

to  check 

radial  intensity  on 

April  20/71 

■ 


. 

H-7 


TABLE  H- 1-6 
RUN  10 


Solution  of  SEPARAN  AP-30  at  a  cone,  of  0,01% 


Q 

=  12.15  USgal/min 

T 

w 

=  0.0656  psf 

T 

=  77,0°F 

* 

U 

=  0.1830  ft/sec 

N'Re 

=  40992 

"V 

=  5.080  ft/sec 

f 

=  0.00282 

<u  > 

x  max 

=  6.12  ft/sec 

n 

y 

inch 

Gx 

ft/sec 

°r 

ft/sec 

Y+ 

U+ 

!x 

V 

1 

5 

0.0056 

1 .686 

-0.003 

7,6 

8.62 

0.105 

0.007 

2 

43 

0.0075 

2,161 

+0.026 

10.2 

11.05 

0.121 

0.010 

3 

62 

0,0115 

2.943 

+0.015 

15.6 

15,05 

0.135 

0.011 

4 

48 

0,0164 

3.402 

+0 » 050 

22,3 

17,40 

0.118 

0,010 

5 

58 

0.0204 

3.580 

+0.044 

27.7 

18,31 

0,124 

0.013 

6 

79 

0.0242 

3.967 

+0,008 

32.9 

20.28 

0,115 

0.012 

7 

79 

0.0309 

4.427 

+0.029 

42.0 

22,64 

0,110 

0.015 

8 

79 

0.0387 

4.724 

+0.028 

52,5 

24.15 

0.109 

0,019 

9 

79 

0.0468 

4.500 

+0.038 

63.5 

23.01 

0,140 

0.022 

10 

79 

0.0623 

4.729 

+0.061 

84.7 

24,18 

0.147 

0.024 

11 

79 

0,0249 

3.933 

+0.020 

30.2 

22.48 

0.117 

0,021 

12 

79 

0,0426 

4.326 

+0 „ 005 

51.8 

24.73 

0,101 

0.023 

13 

79 

0.0648 

4,763 

+0.034 

78.7 

27.23 

0.078 

0.026 

14 

79 

0.1028 

5,009 

+0.043 

123,9 

28,64 

0.080 

0.033 

15 

79 

0.1428 

5.313 

+0.024 

173.5 

30.37 

0,070 

0.028 

16 

79 

0.1826 

5,498 

+0,043 

221.8 

31.43 

0,057 

0.027 

17 

79 

0.2297 

5.637 

+0.077 

279,1 

32,23 

0,060 

0.032 

18 

79 

0.2672 

5.872 

+0 ,038 

325 

33.57 

0.075 

0.016 

19 

55 

0.3012 

5,903 

+0.058 

366 

33.74 

0.066 

0.027 

20 

55 

0.3406 

5,973 

-0.001 

414 

34.15 

0.068 

0.032 

21 

52 

0.4160 

6.018 

+0 , 062 

505 

34.40 

0.061 

0,021 

■ 

. 

Q 

T 

N ' 

Re 

f 

n 

1  18 

2  20 

3  64 

4  79 

5  69 

6  79 

7  79 

8  79 

9  79 

10  79 

11  123 

12  79 

13  110 

14  79 

15  210 

16  94 

17  94 

18  193 

19  210 


H-8 


TABLE  H- 1-7 

RUfJ  11 

water 

=  32.4  USgal/min 

rw 

=  0.753  psf 

=  70.0  °F 

u* 

=  0.624  ft/sec 

=  104299 

<u  > 

X 

=  13.672  ft/sec 

=  0.00442 

<u  > 

x  max 

=  16.40  ft/sec 

y 

i  nch 

ux 

ft/sec 

ur 

ft/sec 

Y+ 

U+ 

I 

r 

0.0159 

10.706 

-0.147 

78.03 

17.17 

0.077 

0.019 

0.0245 

10.897 

+0.218 

120.24 

17.47 

0.050 

0.026 

0.0322 

11 .331 

+0.079 

158.03 

18.18 

0.072 

0.029 

0.0521 

12.301 

+0.068 

255.69 

19.72 

0.063 

0.033 

0.0723 

12.857 

-0.063 

354.82 

20.62 

0.059 

0.033 

0.0906 

13.400 

+0.219 

444.53 

21  .48 

0.062 

0.028 

0.1121 

13.562 

+0.255 

549.90 

21  .75 

0.067 

0.029 

0.1922 

14.696 

-0.005 

943.12 

23.56 

0.064 

0.036 

0.2300 

15.112 

+0.100 

1128.56 

24.23 

0.048 

0.027 

0.3068 

16.013 

-0.225 

1505.86 

25.68 

0.042 

0.036 

0.3490 

15.827 

+0.214 

1712.56 

25.38 

0.047 

0.031 

0.387Q 

16.116 

+0.007 

1899.30 

25.84 

0.044 

0.031 

0.4216 

16.160 

+0.107 

2068.86 

25.91 

0.051 

0 .023 

0.4620 

16.348 

+0 . 1 b6 

2267.32 

26. HI 

0.046 

0.()24 

0.46BD 

16.243 

+0.042 

2296.77 

26.05 

0.039 

0 .030 

0.0906 

13.169 

+0.045 

444.53 

21 .12 

0.068 

0.029 

0.7922 

14.704 

+0.105 

943.10 

23.58 

0.056 

0.033 

0.3870 

15.882 

-0.110 

1899.30 

25.46 

0.048 

0.029 

0.4620 

16.233 

+0.041 

2267.32 

26.03 

0.041 

0.030 

' 

H-9 


TABLE  H- 1-8 
RUN  12 


Solution  of  SEPARAN  AP-30  0.01% 


Q 

=  28.4 

USgal/min 

Tw 

= 

0.2318  psf 

T 

=  74.5  c 

>F 

U  * 

= 

0.3462  ft/sec 

N1 

Re 

=  95817 

<V 

= 

12.128  ft/sec 

f 

=  0.00193 

max 

13.80  ft/sec 

n 

y 

ux 

ur 

Y+ 

u+ 

I 

X 

I 

r 

inch 

ft/sec 

ft/sec 

1 

15 

0.0097 

8.090 

+0.075 

22.3 

23.37 

0.0789 

0.0099 

2 

18 

0.0143 

9.093 

+0.096 

32.8 

26.26 

0.0711 

0.0135 

3 

28 

0.0184 

10.371 

+0.173 

42.2 

29.96 

0.0528 

0.0106 

4 

55 

0.0221 

10.402 

+0.105 

50.8 

30.04 

0.0738 

0.0099 

5 

90 

0.0301 

10.824 

+0.152 

69.1 

31.27 

0.0644 

0.0115 

6 

107 

0.0391 

11.393 

+0.099 

89.7 

32.91 

0.0529 

0.0129 

Q 

=  20.4  USgal/min 

RUN  12B 

Tw 

=  0.2514  psf 

T 

=  74.5  °F 

U* 

=  0.3606  ft/sec 

N  1 

IN  Re 

=  95817 

<u  > 

X 

= 1 2 . 128  ft/sec 

f 

=  0.00179 

<u  > 

x  max 

-  1 3.80  ft/sec 

7 

13 

0.0232 

10.520 

+0.D50 

55.5 

29,17 

0.0514 

0.0162 

8 

28 

0.0311 

10.750 

+0.0192 

74.4 

29.81 

0.0509 

0.0145 

9 

76 

0.0468 

11.335 

+0.022 

112 

31  .43 

0.0558 

0.0137 

10 

100 

0.0665 

11.686 

+0.027 

159 

32.41 

0.0465 

0.0198 

11 

100 

0.1058 

12.373 

+0,054 

253 

34.31 

0,0421 

0,0209 

12 

107 

0.1844 

12.807 

+0,059 

441 

35.52 

0.0394 

0,0181 

13 

97 

0.2630 

13.123 

-0.040 

629 

36.39 

0.0348 

0.0208 

14 

75 

0.3416 

13.462 

-0,047 

817 

37.33 

0.0357 

0.0236 

15 

85 

0.4202 

13.656 

+0.053 

1005 

37,87 

0.0251 

0.0198 

I 


. 

. 

H-1CI 


TABLE  H- 1-9 
RUM  13 

Solution  of  SEPARAN  AP-30  at  a  cone,  of  0.10% 


Q 

=  54.0  USgal/min 

w 

=  0.4682  psf 

T 

=  75.0°F 

* 

U 

=  0.4921 

1  ft/sec 

N' 

Re 

=  37600 

<u  > 

X 

=  21.616  ft/sec 

f 

=  0.00101 

<u  > 

x  max 

=  30.50 

ft/sec 

,,+ 

n 

y 

0 

u 

Y 

U 

I 

I 

X 

r 

X 

r 

inch 

ft/sec 

ft/sec 

1 

72 

0.0101 

5.770 

+0.105 

13.2 

11.72 

0.0491 

0.0039 

2 

57 

0.0140 

7.535 

+0.111 

18.3 

15.31 

0.0349 

0.0030 

3 

63 

0.0180 

9.6D7 

+0.165 

23.5 

19.52 

0,0458 

0,0046 

4 

64 

0.0258 

11.848 

+0.205 

33.8 

24.07 

0.0432 

0.0048 

5 

63 

0.0337 

12.913 

+0.236 

44,1 

26,24 

0.0346 

0.0052 

6 

59 

0.0415 

14.550 

+0.244 

54.4 

29.57 

0.0413 

0.0053 

7 

16 

0.0058 

5.027 

+0.081 

7.73 

10.21 

0.0398 

0.0034 

8 

46 

0.0245 

10.679 

+0.084 

33.5 

21.70 

0.0466 

0.0052 

9 

103 

0.0447 

15.237 

+0.175 

59.3 

30.96 

0.0403 

0.0071 

10 

j 

103 

0.0650 

17.555 

+0.264 

85.5 

35.67 

0.0422 

0.0081 

11 

95 

0.1057 

21.131 

+0.250 

137 

42.94 

0.0438 

0.0086 

12 

95 

0.1452 

22.799 

+0.296 

188 

46.33 

0.0396 

0.0078 

13 

95 

0.1837 

24.714 

+0.298 

240 

50.22 

0.0381 

0.0062 

14 

87 

0.3005 

28.249 

-0.147 

395 

57.40 

0.0427 

0,0121 

15 

87 

0.3796 

29.538 

+0. 388 

498 

60.02 

0.0477 

0,0061 

' 

TABLE  H-I-10 


RUN  14 

Solution  of  SEPARAN  AP-30  at  a  cone.  of  0.10% 


= 

1 1 .274  USgal/min 

= 

77°F 

- 

68450 

= 

0. 00080 

n 

y 

Gr 

inch 

ft/sec 

ft/sec 

105 

0.0119 

+0.249 

12.813 

102 

0.0320 

+0.302 

20.796 

98 

0.0494 

+0.344 

24.623 

103 

0.0711 

+0.310 

27.931 

106 

0.1289 

+0.533 

34.248 

102 

0.2106 

+0.496 

39.365 

12 

0.2494 

+0.421 

40.366 

100 

0.2876 

+0.494 

42.465 

98 

0.3670 

+0.541 

44.83 

100 

0.4839 

0.503 

46.22 

w 

n 

o 

oo 

cn 

o 

psf 

* 

U 

=  0.666 

ft/sec 

<v 

=  33.40 

ft/sec 

<LVmax 

-  49.80 

ft/sec 

Y+ 

U+ 

I 

I 

X 

r 

25.7 

19.25 

0.553 

0.0048 

68.8 

31.24 

0.0529 

0.0054 

106 

36.99 

0.0454 

■0.0073 

153 

41.96 

0.0435 

0.0067 

277 

51.45 

0.0362 

0.0072 

452 

59.14 

0.0351 

0,0077 

536 

60.64 

0.0496 

0.0085 

619 

63.79 

0.0378 

0.0074 

789 

63.75 

0 . 0393 

0.0080 

1041 

69.44 

0.0305 

0.0092 

. 


■ 

H-l  2 


TABLE  H- 1-11 
RUN  15 

water 


Q 

=  93.7  USgal/min 

Tw 

=  0.105  psf 

T 

=  74°F 

* 

U 

=  0.233  ft/sec 

N ' 

V 

=  114083 

‘V 

=  5.00 

ft/sec 

f 

=  0.00435 

<u  >  =6.29 

x  max 

ft/sec 

n 

y 

ux 

ur 

Y+ 

U+ 

!x 

I 

r 

inch 

ft/sec 

ft/sec 

1 

11 

.0106 

2.392 

0.006 

20.40 

10.28 

0.154 

0.009 

2 

58 

.0304 

3.587 

-0.043 

58.51 

15.42 

0.067 

0.038 

3 

78 

.0550 

3.969 

-0.006 

105.9 

17.06 

0.079 

0.046 

4 

95 

.0766 

4.108 

-0.002 

147.4 

17.66 

0.078 

0.038 

5 

95 

.0947 

4.391 

+0.007 

182.3 

18.87 

0.083 

0.051 

6 

95 

.1529 

4.581 

-0.013 

294.3 

19.69 

0.064 

0.049 

7 

95 

.2051 

4.818 

+0.018 

394.8 

20.71 

0.064 

0.041 

8 

95 

.3111 

5.032 

+0.016 

598.8 

21.63 

0.064 

0.038 

9 

95 

.4099 

5.320 

+0.045 

789.0 

22.87 

0.060 

0.033 

10 

45 

.6071 

5.641 

-0.014 

1168 

24.24 

0.072 

0.052 

n 

42 

.8021 

5.790 

+0.008 

1543 

24.89 

0.051 

0.031 

12 

95 

.702 

5.977 

+0.004 

1351 

25.69 

0.053 

0.036 

13 

55 

1.00 

6.285 

-0.024 

1925 

27.01 

0.056 

0.030 

14 

26 

1.20 

6.240 

-0.092 

2310 

26.82 

0.058 

0.028 

15 

38 

1 .40 

6.290 

+0.009 

2695 

27.03 

0.053 

0.022 

. 

H-l  3 


TABLE  H- 1-1 2 
RUN  16 

Solution  of  SEPARAN  AP-30  at  a  cone,  of  0.01% 


Q 

T 

N, 

f 


1 

2 

3 

4 

5 

6 

7 

8 
9 


— 

95  USgal/min 

Tw 

=  0.065  psf 

= 

74°F 

* 

U 

=  0.183  ft/sec 

= 

97774 

A 

xC| 

V 

=  5.07 

ft/sec 

— 

0.00258 

<u  >  =  5.87 

x  max 

ft/sec 

n 

y 

inch 

ux 

ft/sec 

Ur 

ft/sec 

Y+ 

U+ 

!x 

lr 

89 

.0252 

3.458 

0.032 

31.2 

18.88 

0.112 

0.012 

95 

.0528 

3.695 

0.023 

85.4 

20.17 

0.089 

0.024 

95 

.1017 

4.368 

0.048 

126 

23.84 

0.067 

0.027 

95 

.2014 

4.806 

0.034 

249 

26.23 

0.057 

0.027 

95 

.300 

4.989 

0.050 

372 

27.23 

0.048 

0.028 

95 

.505 

5.340 

0.068 

626 

29.15 

0.046 

0.033 

95 

.702 

5.531 

0.069 

870 

30.19 

0.043 

0.023 

75 

.905 

5.715 

-0.053 

1121 

31 .19 

0.033 

0.024 

70 

1.101 

5.798 

0.060 

1365 

31 .65 

0.043 

0.026 

90 

1.30 

5.844 

0.066 

1611 

31 .90 

0.033 

0.022 

10 


✓ 
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TABLE  H-I-13 
RUN  17 

Solution  of  SEPARAN  AP-30  at  a  cone,  of  0.10% 


Q 

=  163  USgal/min 

Tw 

=  0.106  psf 

T 

=  74°F 

* 

U 

=  0.234  ft/sec 

N  1 
NRe 

=  21905 

"V 

=  8.77 

ft/sec 

f 

=  0.00146 

<uxW1K38 

ft/sec 

n 

y 

Ux 

ur 

Y+ 

U+ 

*x 

I 

r 

inch 

ft/sec 

ft/sec 

1 

22 

0.0204 

1.671 

0.049 

6.27 

7.12 

0.041 

0.0064 

2 

60 

0.0514 

3.117 

0.038 

15.4 

13.29 

0.064 

0.0069 

3 

91 

0.1018 

5.069 

0.074 

30.4 

21 .61 

0.065 

0.0097 

4 

95 

0.1565 

6.714 

0.100 

47.0 

28.62 

0.074 

0.0090 

5 

87 

0.2043 

7.364 

0.087 

61.3 

31  .39 

0.072 

0.0098 

6 

95 

0.3023 

8.863 

0.088 

90.6 

37.78 

0.059 

0.0095 

7 

79 

0.4004 

9.7229 

0.124 

119 

41.45 

0.047 

0.0109 

8 

84 

0.6005 

10.524 

0.120 

180 

44.86 

0.044 

0.0108 

9 

30 

0.7971 

10.989 

0.029 

239 

46.85 

0.041 

0.114 

10 

33 

0.999 

11.409 

0.046 

300 

48.64 

0.035 

0.0084 

11 

35 

1.1985 

11.252 

0.113 

359 

47.97 

0.025 

0.0110 

12 

35 

1.298 

11.360 

0.086 

389 

48.43 

0.034 

0.0103 

' 

• 

APPENDIX  I 


1-1 


PRESSURE  DROP  MEASUREMENTS 


Pressure  drops  over  a  range  of  flow  rates  were  obtained 
for  each  of  the  experimental  fluids  and  are  tabulated  in  Tables  I-I. 
For  the  Newtonian  fluids  the  data  serve  to  validate  the  apparatus 
while  for  the  drag  reducing  solutions  the  data  characterize  the 
drag  reduction  over  a  range  of  flow  rates  of  interest.  The 
tables  representing  pressure  drop  measurements  during  a  photographic 
run  are  specified  so  that  degradation  effects  may  be  following. 


1-2 


TABLE  I-I-l 

FRICTION  COEFFICIENT  -  GENERALIZED  REYNOLDS  NUMBER 

Water 


Pipe  Diameter  =  0.833  ft 


Length 

Between 

Pressure  Taps  = 

11  ft 

Flow 

T 

Manometer 

Friction 

General i zed 

Friction 

Readi ngs 

loss 

Reynolds 

Factor 

Gal /Min 

°F 

Inches 

psf 

Hg 

43.2 

84.5 

10.38 

668.0 

164200 

0.00423 

42.0 

72.0 

9.39 

612.0 

138000 

0.00410 

21 .2 

68.0 

2.75 

181.2 

66000 

0.00469 

16.3 

68.0 

1.84 

121 .2 

50700 

0.00529 

18.7 

74.4 

2.30 

149.5 

64000 

0.00459 

19.7 

74.7 

2.50 

162.5 

67400 

0.00484 

9.5 

76.0 

0.81 

52.7 

17200 

0.00675 

' 

; 

1-3 


TABLE  I- 1-2 

FRICTION  COEFFICIENT  -  GENERALIZED  REYNOLDS  NUMBER 


Water 

Pipe  Diameter  =  .2300  ft 


Temperature  = 

70°F 

Flow 

Manometer 

Length 

General i zed 

Fri cti on 

Readi ngs 

(Taps) 

Reynolds 

Factor 

Gal /Min 

Inches 

Ft 

ODB 

113.2 

39  „  60 

23.60 

131776 

0.0043 

113.2 

29.45 

17.50 

132014 

0.0043 

113.2 

19.12 

11,69 

132150 

0.0042 

96.5 

29 . 65 

23.60 

112336 

0.0044 

96.5 

21 .90 

17.50 

112538 

0,0044 

96.5 

14.42 

11.69 

112654 

0.0043 

83.5 

22.65 

23.60 

97202 

0,0045 

83.5 

16.81 

17,50 

97378 

0.0045 

83.5 

11.10 

11,69 

97478 

0.0044 

70.5 

16.86 

23.60 

82069 

0,0047 

70.5 

12.55 

17.50 

82217 

0.0047 

70.5 

8.30 

11,69 

82302 

0.0046 

61 .2 

12.73 

23,60 

71243 

0,0047 

61.2 

9.52 

17.50 

71371 

0,0047 

61 .2 

6.32 

11,69 

71445 

0.0047 

51 .2 

9.70 

23.60 

59602 

0.0052 

51 .2 

7.25 

17,50 

59709 

0,0052 

51 .2 

4.89 

11.69 

59771 

0.0052 

43.3 

7.10 

23.60 

50406 

0.0053 

43.3 

5.27 

17.50 

50496 

0,0053 

34.0 

4.62 

23.60 

39579 

0,0056 

34.0 

3.54 

17.50 

39651 

0.0057 

Water 

Pipe  Diameter 

=  .2300  ft 

Temperature  74 

°F 

Flow 

Manometer 

Length 

General i zed 

Friction 

Readi ngs 

(Taps) 

Reynolds 

Factor 

Gal /Min 

Inches 

Ft 

Hg 

92.0 

0.66 

23,60 

112645 

0.0045 

92.0 

0.50 

17,50 

112847 

0.0045 

126.0 

1.16 

23.60 

154274 

0.0042 

126.0 

0.87 

17,50 

154552 

0.0042 

142.5 

1  .45 

23,60 

174477 

0.0041 

142.5 

1.10 

17.50 

174791 

0.0041 

■ 

TABLE— I— I—  3  FRICTION  COEFFICIENT  -  GENERALIZED  REYNOLDS  NUMBER 


1-4 


z 

o 


#— « 

rv 

vf 

o 

00 

0D 

m 

r\j 

m 

« — 4 

a) 

nj 

CO 

in 

m 

00 

go 

O' 

gf 

in 

vT 

oo 
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